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ABSTRACT 


The  general  equations  governing  the  angular  motion  of  a  rigid 
satellite  are  investigated  for  the  case  where  both  gravitational  and 
aerodynamic  torques  are  included,  and  ^era  the  atmosphere  is  assumed  to 
be  spherically  symmetric  but  rotating.  For  a  non-spinning,  axisymmetric 
satellite,  these  equations  show  the  vehicle  to  be  characterized  by  three 
dimensionless  aerodynamic  parameters  and  one  inertia  parameter.  A 
Fourier  series  expansion  of  the  coefficients  of  these  equations  is  found 
to  yield  equations  of  the  non-homogeneous  Hill  type. 

The  complete  equations  of  motion  are  solved  for  the  case  of  a 
circular  orbit.  It  is  shown  that  atmospheric  rotation  induces  a  steady- 
state  yawing  motion  irtiich  prevents  yaw  equilibrium  from  being  achieved 
for  other  than  equatorial  orbits.  The  effects  of  configuration  and 
orbit  parameters  are  examined  in  detail.  It  is  found  that  gravity 
gradients  will  not  generally  dominate  tuitil  altitudes  in  excess  of  300- 
UOO  miles  are  achieved,  while  at  lower  altitudes,  even  small  amounts  of 
aerodynamic  stability  suffice  to  override  the  gravitational  torques. 

The  stability  of  the  pitching  motion  on  an  eccentric  orbit  is 
analyzed.  It  is  found  that  an  infinity  of  unstable  regions  may  be 
encountered  in  which  small  deviations  from  a  purely  periodic  motion  will 
amplify.  Computer  results  are  presented  to  Illustrate  the  catastrophic 
nature  of  these  unstable  regions  vrtiich  may  appear  in  a  situation  giving 
the  outward  appearance  of  stability.  An  approximate  solution  for  the 
pitching  motion  on  an  elliptical  orbit  is  also  obtained  using  the  method 
of  variation  of  parameters.  The  results  are  compared  with  a  numerical 
integration  of  the  complete  governing  equation. 

It  is  generally  concluded,  however,  that  the  aerocfynamic  and 
inertial  design  of  a  satellite  may  be  acconplished  so  as  to  minimize  or 
altogether  eliminate  the  unstable  regions.  In  particular,  the  design  of 
the  vehicle  for  maximum  aeroc^andc  stability  appears  to  have  proodse 
for  achieving  wholly  passive  attitude  stabilization  at  low  altitudes. 
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INTRODUCTION 


The  problem  of  achieving  precise  attitude  stabilization  of  artifi¬ 
cial  satellites  has,  of  late,  received  increasing  attention,  principally 
because  adec[uate  means  of  orientation  and  control  are  necessary  prerequi¬ 
sites  to  the  serious  scientific  and/or  military  use  of  such  vehicles. 

The  general  subject  may  be  divided  into  two  rather  distinct  but  inter¬ 
related  areas:  (l)  the  determination  of  the  external  torques  which  act 
on  the  satellite  together  with  the  passive  response  of  the  vehicle  to 
various  combinations  of  these  torques,  and  (2)  the  synthesis  of  active 
attitude  sensing  and  control  systems  idiich  utilize  these  torques  to 
achieve  and  hold  some  desired  vehicle  orientation. 

The  general  framework  within  which  this  problem  may  be  treated  has 
been  discussed  by  numerous  authors,  and  reasonably  comprehensive  surveys 
of  the  field  have  been  published  by  Roberson  (Ref.  3)  and  Frye  and 
Stearns  (Ref.  li) .  Since  the  synthesis  of  attitude  control  systems  can¬ 
not,  in  general,  be  accomplished  without  detailed  knowledge  of  the 
vehicle's  response  to  the  applied  external  torques,  this  latter  subject 
has  received  the  most  extensive  treatment.  The  majority  of  the  published 
work  in  this  area  has,  however,  been  confined  to  examining  the  effect  of 
only  a  single  external  torque  and/or  the  assumption  of  a  circular  orbit. 

For  example,  the  effects  of  gravity  gradient  torques  alone  on  the 
librations  of  a  dumbbell-shaped  satellite  have  been  studied  by  Schindler 
(Ref.  3)  and  Klemperer  and  Baker  (Ref.  6)  for  the  circular  orbit  case, 
while  Stocker  and  Vachino  (Ref.  7)  and  Baker  (Ref.  8)  have  discussed  the 
motion  on  an  eccentric  orbit.  Roberson  (Refs.  9,  10)  has  qualitatively 
discussed  this  entire  problem  with  considerable  generality. 

The  influence  of  aeroctjmamic  torques  on  a  vehicle  in  a  circular 
orbit  has  been  examined  by  Wall  (Ref.  11)  and  the  combined  effect  of 
aerodynamic  and  gravitational  torques  has  been  covered  briefly  by  DeBra 
(Ref.  12)  and  DeBra  and  Stearns  (Ref.  13),  again  for  the  circular  orbit 
case.  Beletskiy  (Ref.  Ih)  has  also  considered  this  problem  in  a  some¬ 
what  more  general  form.  Along  similar  lines,  the  solar  torques  have  been 
discussed  by  Sohn  (Ref.  15),  Frye  and  Stearns  (Ref.  U)  euid,  more  recently, 
by  Newton  (Ref.  16). 

Almost  without  exception,  none  of  the  papers  have  attempted  to  evolve 
rational  design  criteria  or  to  provide  detailed  analytical  results  for  the 
vehicle's  angular  motion  when  combinations  of  several  torques  act.  In 
particuljur,  the  combined  Influence  of  aerodynamic  and  gra'^tational 
torques  has  apparently  not  been  carefully  considered  in  spite  of  the  fact 
that  these  may  be  readily  shown  to  be  the  dominant  external  torques  for 
orbits  below  about  hOO  miles  altitude.  Although  DeBra  (Ref.  12)  draws 
some  preliminary  conclusions  about  the  interaction  of  these  two  torques. 
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he  has  considered  only  the  equilibrium  case  and  has  greatly  iver simplified 
the  aerc-^— •>ric  treatment.  Moreover,  no  analytical  consideration  of 
ellipt' ,  bits  has  been  made. 

Since  the  satellite's  external  configuration  determines  the  aero¬ 
dynamic  torques  in  a  manner  analogous  to  the  way  the  mass  distribution 
determines  the  gravitational  torques,  it  is  clear  that  an  understanding  of 
both  of  these  torques,  and  .their  mutual  interaction,  is  essential  for  a 
rational  design  to  be  made . 

Accordingly,  this  work  seeks  to  extend  the  theory  of  satellite  atti¬ 
tude  control  by  examining  the  detailed  aerodynamic  influence  on  the 
angular  motion  including  interactions  with  the  gravitational  torques  and 
considering  a  general  elliptical  orbit  of  small  eccentricity.  The  effect 
of  atmospheric  rotation  is  also  considered  to  the  extent  to  which  it 
affects  the  vehicle's  angular  motion.  Torques  due  to  radiation  pressure, 
electromagnetic  effects  and  internal  moving  parts  are  considered  briefly 
but  are  concluded  to  be  negligible  over  the  altitude  range  for  which  aero¬ 
dynamic  effects  are  significant.  The  motion  of  the  vehicle's  mass  center 
(which  is  only  weakly  coupled  to  the  angular  motion  through  the  gravita¬ 
tional  gradient  and  the  drag)  is  assumed  to  be  independent  and  specified. 
As  shown  by  Moran  (Ref.  29),  the  uncoupling  thus  performed  introduces 
essentially  no  error  into  the  final  results  for  the  librations,  at  least 
for  small  eccentricity  orbits. 
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1.  EQUA'i’IONS  OF  MOTION  AND  COORDINATE  SYSTEMS 


1.1  Complete  EQuations  of  Motion 

For  a  rigid  body  of  constant  mass  possessing  principal  axes  of 
inertia  ’*■  j  ^  equations  of  angular  motion  about  the  mass 

center  are  the  well-known  Euler  Equations  of  a  rigid  rotator: 

L.  =  Tl.  -//»  (VI,) 

where  'f  ,  t  ,  tl  are  components  of  the  instantaneous  angular  velocity 
of  the  body  (relative  to  an  inertial  reference  frame)  along  the  principal 
body-fixed  axes  and  L*.  ,  ^“4  ’  '’'^®  corresponding  components  of 

the  instantaneous  applied  external  torque.  The  products  of  inertia  are 
absent,  of  course,  due  to  the  choice  of  pi-incipal  axes.  The  applied 
external  torques  will  be  examined  in  detail  in  Section  2.  The  motion  of 
the  body's  mass  center  is  assvimed  to  be  known  and  independent  of  the 
angular  motions  about  the  mass  center. 

The  problem  is  to  choose  an  orientation  scheme  such  that  the  equa¬ 
tions  of  motion,  Eqs.  (1),  assume  their  simplest  form.  Since  the  method 
of  orientation  determines  the  simplicity  or  complexity  of  the  angular 
velocity  vector,  judicious  choice  of  coordinates  is  of  considerable  value. 
Roberson  (Ref.  9)  has  discussed  this  problem  qualitatively  at  some  length. 
Although  only  three  independent  coordinates  are  required,  it  is  desirable 
to  introduce  more  than  this  number  in  order  to  retain  a  clear  physical 
significance  for  each  coordinate .  The  inertial  reference  frame  chosen 
here  is  oriented  with  respect  to  the  planet  as  shown  in  Fig.  1.  The 
"Y  -axis  lies  along  the  planet's  axis  of  rotation  and  points  south,  while 
X  and  2^  lie  in  the  equatorial  plane  with  —  %  directed  toward  the  most 
recent  vernal  equinox  in  the  customary  astronomical  fashion^. 


lit  is  assumed  here  that  the  only  significant  motion  of  the  earth  which 
precludes  its  use  as  an  inertial  reference  frame  is  its  daily  rotation. 
Thus,  the  inertial  coordinates  defined  above  differ  from  geocentric 
earth-fixed  coordinates  by  a  uniform  rotation  of  the  latter  about  their 
common  • 


Follovdng  the  astronomical  convention  (Ref.  17),  three  Euler  angles 
are  adopted  to  orient  the  instantaneous  orbit  and  orbit  plane  relative  to 
the  inertial  reference  frame.  These  are  the  nodal  longitude  (or  right 
ascension)  ,  the  orbital  inclination  and  the  argument  of  perigee, 
tO  .  Their  general  relationship  is  shown  in  Fig.  1,  where  it  is  noted 
that,  owing  to  the  choice  of  inertial  coordinates,  the  rotations  ,  i 

and  CO  are  left-handed.  The  radial  distance  from  the  planet's  mass 
center  to  the  vehicle's  mass  center  is  denoted  by  r  . 


Figure  1.  Orientation  of  the  Orbit  Plane  in  Space 


The  result  of  performing  these  rotations  is  a  right-handed,  orthogonal 
coordinate  triad  with  origin  at  the  primary  focus  of  the  instantaneous 
ellipse  with  pointing  toward  perigee  and  X’p  pointing  in  the 

direction  of  vehicle  motion  through  perigee  (Fig.  1).  In  the  absence  of 
all  perturbations,  the  vehicle's  motion  is  given  by  the  well-known 
Keplerian  laws  of  motion  and  the  angles  ^  ,  C  and  u?  are  constants. 
When  perturbative  influences  such  as  those  due  to  planetary  oblateness, 
atmospheric  drag,  etc.,  are  considered  these  angles  change  with  time, 
having  secular  and/or  periodic  variations,  but  the  rates  are  generally 
small  (Ref.  18).  These  angular  rates  will  be  considered  further  shortly 
and  will  be  shown  to  have  only  a  second-order  effect  on  the  vehicle's 
motion. 


■1  • 

■‘‘It  is  conventional  to  denote  the  orbit  inclination  by^  ;  however  i  is 

used  here  to  avoid  later  confusion. 


The  location  of  the  vehicle  on  the  orbit  path  is  most  conveniently 
specified  by  the  tirue  anomaly,  ir  ,  which  is  measured  in  the  orbit  plane 
from  perigee  and  is  positive  in  the  direction  of  motion.  The  angles 
^2-  ,  L  and  aJ-t-xTaa  define  a  set  of  orbital  coordinates  which  move 
with  the  vehicle,  and  are  con^letely  determined  by  the  translational 
motion  of  the  vehicle's  mass  center.  The  transformation  from  the  inertial 
coordinates  to  tne  orbital  coordinates  is  given  in  matrix  notation  by 
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The  abbreviations  C  =  cosine  and  S  *  sine  have  been  used  here  for 
brevity.  The  appropriate  geometry  is  shown  in  Fig.  2.  The  axis, 

which  is  always  directed  inward  along  the  radius  vector,  represents  a 
sort  of  "dovm"  direction,  while  x,  ,  which  points  in  the  direction  of 
vehicle  motion,  indicates  the  "forward"  direction. 
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Figure  2.  Orbital  Coordinate  System 


Inatx;.;uch  as  the  orbital  coordinates  are  conpletely  defined  by  the 
:r.otion  of  the  vehicle's  mass  center,  three  independent  angles  are  yet 
required  to  orient  the  vehicle.  The  choice  is  again  dictated  by  the 
resulting  simplicity  of  the  e>^aations  of  motion.  Roberson  (Ref.  9)  and 
DeBra  (Ref.  12)  have  pointed  out  that  >^en  gravitational  torques  are  to 
be  considered,  the  first  of  these  three  remaining  angles  should  be  formed 
by  a  rotation  about  >  since  the  gravity  gradient  (for  a  spherical 
central  body)  occurs  along  this  axis.  This  point  is  also  demonstrated  in 
Section  2.2.  On  the  other  hand,  for  the  current  work  where  aerocjynamic 
torques  are  of  interest,  it  is  desirable  to  keep  the  first  angle  in  the 
nominal  pitch  plane  and  so  the  first  rotation  should  occur  about  Yo  • 
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Roberson,  in  his  ext/ensive  treatment  of  this  subject  (Refs.  3,  9,  10) 
has  shora  that  for  small  angles,  the  equations  of  motion  uncouple  such 
that  the  order  in  vdii.ch  the  rotations  are  performed  is  immaterial.  However, 
since  it  may  be  desir  ’  .le  to  examine  larger  deviation  angles  in  pitch^,  the 
first  rotation  about  Yo  will  be  employed.  The  angles  which  shall  be  used 
here  are  denoted  by  ©  ,  y-  ,  ^  and  are  defined  as  shovm  in  Fig.  3.  Thus 
the  transformation  from  the  orbital  coordinates  to  the  body  axes  is 


where 
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^For  a  non-spinning  vehicle,  Roberson  (Ref.  3)  has  shown  the  pitching 
motion  to  be  independent  of  the  rolling  and  yawing  motions,  vrtiich  are 
inherently  coupled.  This  point  is  also  discussed  in  Section  l.lj. 


Figure  3-  Vehicle  Orientation  Angles 


The  complete  orientation  transformation  from  the  inertial  coordinates 
to  the  principal  body  axes  is  therefore  composed  of  six  angles  as  follows: 


~  X  ' 

"x 

Y 

(6) 
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Here  &  ,  and  are  the  only  unknowns,  since  the  center  of  mass 

motion  is  assumed  to  be  specified.  The  principal  body  axes  defined  above 
are  shown  in  relation  to  the  vehicle  in  Fig.  1^. 


Figure  U.  Vehicle  Body  -  Fixed  Coordinate  System 


The  angular  velocity  components  in  the  body  axes  are  readily  found  from 
the  transformation  equation,  Eq.  (6),  in  the  following  general  form; 
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where  s  is  a  possible  spin  velocity  about  the  <  axis.  It  is  noted, 
however,  that  it  is  convenient  to  define 


s  ©  -  U 
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since 
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Expanding  Eq.  (8)  there  is  obtained, 

T  ~  i>  +S  +£  Svntrjc^Y'  +^[s«n-^C^6  +  ij 

J  =  -  y-  Sm  ^  C^<f>  Co-.y  -  Co.^^  Sin  y  Sinl^^ 

+<Si[sin^Sir«t;fSint  -  C«>-yCcM»  C  +-  Co^j^SioV  j 

=  — y>  —  t|  Sin^  Co-»y  — Co-<f  f  Svn^  sin  ^ 

+  ^  ^C.«^^Scn|^  Sin  d  +  Sin^  Sin^  Sini^  Sin 

which  are  the  complete  angular  velocity  components  in  the  body  axes. 
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1.2  Effects  of  Orbit  Perturbations 


Before  proceeding  with  the  expansion  of  the  various  torque  components, 
it  is  desirable  to  simplify  the  equations  of  motion  somewhat,  by  simplify¬ 
ing  the  angular  velocity  components  T  ,  ^  ,  h,  •  For  order  of  magnitude 
purposes,  all  angular  velocities  may  be  compared  with  the  angular  velocity 
of  the  orbital  coordinates,  v-  . 

•  •  • 

As  mentioned  previously,  the  angular  velocities  ef2  ,  ^  and  to 
arise  because  of  perturbations  from  the  otherwise  inverse-square  (Kepler- 
ian)  motion.  Examples  are  anomalies  in  the  gravitational  field  (such  as 
oblateness)  and  atmospheric  drag.  These  perturbations  affect  not  only 
the  angles  <SZ.  ,  I  and  co  but  also  the  orbital  period,  the  perigee 
distance  and  other  parameters  which  characterize  motion  on  the  Keplerian 
ellipse  (Ref.  17).  The  perturbations  are  of  two  types,  secular  perturba¬ 
tions,  which  are  monotonous  functions  of  time,  and  periodic  perturbations, 
which  vary  harmonically  with  (say)  the  true  anomaly,  .  For  motion  near 
an  oblate  primary  body,  the  gravitational  potential,  $  ,  may  be  expanded 
as  a  series  in  spherical  harmonics,  the  coefficients  of  which  are  functions 
of  the  radial  distance  r  , 

i  =-=f{  '  ^  ■"■■■] 

where  is  the  product  of  the  universal  gravitational  constant  and  the 
mass  of  the  primary  body,  9c  is  the  mean  equatorial  radius  of  the 
primary  body,  2^^  the  geocentric  latitude  measured  from  the  equator  and 
ij  is  a  small  constant.  Using  the  standard  methods  of  variation  of 
parameters  (Ref.  17),  Krause  (Ref.  18)  finds  that  when  J  is  not  negli¬ 
gible,  the  leading  perturbation  terms  for  and  are  secular  while 
I  has  no  first  order  secular  component,  but  experiences  an  appreciable 
periodic  perturbation.  These  results  may  be  written  in  the  form: 


k=-  - 


m 

CO  = 


where  ^  is  the  orbital  eccentricity,  a  the  mean  distance  (or  semi¬ 
major  axis)  of  the  elliptical  orbit  and  -n  the  Keplerian  mean  motion 
(or  mean  angular  velocity  of  the  motion).  For  nearly  circular  orbits 
(Ref.  18), 


L 


J  -  \ 

+  6c«-er\rj  "Z-U.  Sun  4 


(10) 
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where  p  is  the  orbital  period,  and  hence  for  order  of  magnitude  purposes 
it  is  approximately  true  that 


"(i)®«(i)=  ‘'(i) « j- 

For  Earth,  J  =  I.&  >^io  *  and  hence  the  above  terms  are  negligibly  small 
con^jared  to  unity.  Beletskiy  (Ref.  lU)  arrives  at  essentially  this  same 
result  by  means  of  substantially  similar  arguments. 

The  perturbations  due  to  atmospheric  drag  are  not  so  readily  dis¬ 
posed  of.  Krause  (Ref.  l8)  has  shown  that  when  the  planetary  atmosphere 
does  not  rotate,  drag  produces  no  change  at  all  in  the  angles  JZ.  and 
£  and  so  the  orbit  remains  on  a  space-fixed  plane.  Under  these  circum¬ 
stances  tc>  possesses  a  small  periodic  perturbation  vrtiich  may  be  taken 
into  account  through  the  variable  1^  .  More  recently,  Sterne  (Ref.  19) 
has  considered  the  problem  in  great  generality,  accounting  for  atmospheric 
rotation  as  well  as  the  "flattening"  of  the  atmosphere  due  to  planetary 
oblateness.  In  the  presence  of  atmospheric  rotation,  Sterne's  results  for 
the  secular  pertvirbations  of  ,  C  and  u?  may  be  written  approxi¬ 
mately  (for  small  eccentribity  orbits)  as 


4  = 


CO 

V- 


.a 


Co-»  c 


where  b  is  the  ballistic  coefficient  of  the  bo4y,  cOg  is  the  Earth's 
angular  velocity  of  daily  rotation,  the  atmospheric  density  at  the 

semi-latus  rectum  of  the  orbit  and  1„(,X)  is  the  modified  Bessel  function 
of  the  first  kind  and  order  vt  of  argument  X  .  Here  X  is  an  atmos¬ 
phere  parameter  which  is  defined  in  Section  3.2.  It  will  be  shown  later 
that  for  Earth,  (u3b/-u-)  <.  o(to“*)'trtiile  aR<  £■  O  (icT^)  and  since  b 
is  generally  at  most  of  order  10,  then  the  secular  contributions  to  di  , 

£  and  (O  due  to  atmospheric  drag  on  the  satellite  are  also  negligible 
in  comparison  to  v  . 
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The  resulting  anj^lar  velocity  components  are  found  from  Eqs.  (9)  as, 

•  , 

?  •*  ^  +  S  ~  SinV- 

^  —  -p  <f>  +  tj  Co-t-p  (U) 

^  ss  —p  —  ty  Sin^  C.o-a'j 

where  the  equations  of  the  angular  motion  are  given  by  Eq.  (1)  and  are 
complete  once  the  torque  components  L*  ,  L.^  and  are  specified. 

1.3  Introduction  of  Gravity  Gradient  Torques 

The  determination  of  the  external  torques  acting  on  the  satellite 
will  be  performed  in  Section  2.  However,  it  is  noted  that  one  primary 
torque  known  to  act  on  any  satellite  is  that  due  to  the  gravity  gradient 
across  the  vehicle.  Since  this  torque  arises  from  the  vehicle's  mass 
distribution,  it  bears  a  close  association  with  the  body's  inertia  charac¬ 
teristics  and  is  therefore  best  included  in  the  basic  equations  of  motion. 

From  Eq.  (A-15),  Appendix  A,  the  gravitational  torques  for  an 
inverse-square  central  force  field  may  be  given  to  a  first  order  in  the 
ratio  -^/r  (where  Jl  is  the  maximum  body  dimension)  by 
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where  ^  is  the  product  of  the  earth's  mass  and  the  universal  gravita¬ 
tional  constant  as  before  and  r»/ r  ,  r^/r  and  'T^./r  are  direction 
cosines  between  the  radius  vector  and  the  *  ,  ,  5.  bocfy-fixed  axes. 

These  direction  cosines  may  be  expanded  in  terms  of  the  angles  ©  ,  V  , 
p  from  Eqs.  (U)  and  (5)  as 
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and  so. 


■M  «• 

4* 

[s  2  ^  (s*^s*©  -  c^)  +  -  sV)  s  ](I^  - 1 J 

4^ 

“  z  r» 

[s^sz-^  s»©  +  c^cVsz© 

[-  S?5  c^ts^® 

(Hi) 


The  complete  equations  of  angular  motion  may  now  be  obtained  by 
substituting  Eqs.  (lU)  into  Eqs.  (1).  Because  of  the  aforementioned 
relationship  between  the  gravitational  torques  and  the  angular  inertia, 
it  is  convenient  to  group  these  terms  together,  viz: 


f  L  ^-•|^[s^^CsVs•9-c■‘©) 

=  i  siiJ  5‘e + I  (i^-  i^) 

where  ,  Le^,  denote  all  other  torques  acting  on  the 

vehicle  except  those  due  to  the  gravity  gradients. 


Eqs.  (15)  constitute  the  complete  equations  of  angular  motion  of  a 
"small"  satellite  vehicle  moving  in  the  vicinity  of  an  essentially 
spherical  primary  body.  The  equations  are  written  in  principal  body  auces, 
but  apart  from  these  restrictions,  no  other  significant  assumptions  have 
been  made. 


l.U  Equations  of  Motion  for  Small  Angles 

Now  suppose  that  j>  and  are  small,  but  &  need  not  be  small. 
Then  squares  and  products  of  ^  and  -jA  are  negligible  but  sines  and 
cosines  of  &  must  be  retained.  Eqs.  (l5)  became,  under  these  conditions, 

=  f  I*  +3^['}^s.fl©c^© 

=  i  1^  +{t/1  +3^S.nec^©](!^-!p 

+3^[^s»n©a«©  --54sin'©]|(I^-Ix) 


(16) 


where  the  jmgular  velocity  components  from  Eqs.  (11)  simplify  to  read, 


t  =  4  +-S  - 

I  S  (17) 

It  =  -4  • 

•  • 

To  the  same  order  of  approximation,  if  the  rates  V*  and  (f>  are  small  so 
that  squares  and  products  of  these  terms  are  negligible  also, 

in  S--5^ti 
yti  =-s(;^^+-^) 

ri  =  +Cs+^)»i 

and  (if  s  =  constant), 

•f  =  4  -•V'fi 

/  f[  (18) 

.i  S 

whereupon,  Eqs.  (16)  take  the  form, 

n  1^  -  ^sC^fj+T^)  -3^sm&c^6>j  (19) 

These  equations  describe  the  angular  motion  of  a  satellite  vehicle  in 
the  presence  of  a  non-negligible  spin  velocity,  allowing  &  to  be 
arbitrarily  large.  For  most  applications  of  interest,  however,  ^  is 
also  small  and  in  this  case  the  equations  become 
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'-R*  I'J  ’■^x)  -  ^('i 

nT.j,-(s(W+lf-)-3=^e](i,-lj)  (20) 

These  equations  are  essentially  identical  to  those  presented  by  Fiye  aind 
Steams  (Ref.  U)  and  Roberson  (Ref.  9)>  although  these  authors  do  not 
consider  the  possibility  of  a  spinning  vehicle.  When  the  spin,  s  ,  is 
zero,  the  second  of  Eqs.  (20)  is  uncoupled  from  the  first  and  third  and, 
provided  that  is  independent  of  ^  and  'p'  ,  the  &  variation 

(and  hence  the  rj  Variation)  may  be  explicitly  determined.  The  first  and 
third  of  these  equations  then  contain  only  the  terms  -p-  and  p  as 
coupled  \mknowns  and  may  be  solved  by  iteration.  Note  further  that  since 
the  time  variation  of  rj  may  be  found  from  the  second  of  Eqs.  (20),  the 
above  system  of  equations  is  linearl,  in  spite  of  the  obvious  complexity. 

For  simplicity,  only  the  case  of  a  non-spinning,  S3nnmetrical  vehicle 
will  be  considered  in  detail  herej  that  is,  one  for  which  the  x-axis  is 
an  axis  of  symmetry; 


and  for  which  S  =o.  The  equations  of  angular  motion  become,  in  this 
case, 

=  'll  +3^®(ix-i)  (21) 

=  -/nlx'+T'-iHix-i)  . 

It  is  to  be  noted  from  Eqs.  (20)  that  for  small  angles,  the  gravita- 
ticmal  torques  have  no  direct  effect  on  the  yawing  motion.  For  a  body  of 
revolution  (1^=»X^  =  X^  >  these  torques  also  vanish  from  the  rolling 


L 


and 


must  also  be  linear. 
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moilon  equation  as  shown  by  Eqs.  (21).  The  gravitational  torques  do, 
however,  influence  these  notions  through  their  coupling  with  the  pitching 
motion . 

Finally,  it  is  appropriate  to  employ  the  definition  of  ,  Eq.  (7), 
and  note  that  for  small  angular  rates. 


S  -  a 

^  =  —  u,  5$ 


Hence  Eqs.  (21)  become,  after  some  algebraic  rearrangement, 


•  •  • 

<f>  =  V 


G  +  3  ^  M  e  -  /i  ^  w- 


(22) 


where  M  is  an  inertia  parameter,  defined  by 


M  H 


(23) 


Eqs.  (22)  are  the  equations  of  angular’ motion  to  be  considered  in 
detail.  Inasmuch  as  the  orbital  motion  of  the  vehicle  is  assumed  to  be 
known,  U.  and  all  of  its  derivatives  are  known  and  Eqs.  (22)  ms^  be 
solved  for  the  three  unknowns  ^  as  functions  of  time.  It 

is  well  to  re-emphasize  that  the  rolling  and  yawing  motions  are  inherently 
coupled,  principally  because  u  is  not  negligible  (see  Eqs.  22)j  that  is, 
because  the  orbital  coordinates  are  rotating  with  angular  velocity  tt 
relative  to  the  inertial  axes.  Roberson  (Ref.  3)  has  also  made  this 
observation  previously. 

Before  seeking  a  solution  to  Eqs.  (22),  the  expressions  for  the 
remaining  moments  acting  on  the  vehicle  must  be  obtained.  This  develop¬ 
ment  is  performed  in  the  next  section. 


2 .  EXTERNAL  TORQUES  ON  THE  SATELLITE 


2.1  Introduction 


One  can  readily  conceive  a  number  of  torques  viiich  may  act  on  an 
orbital  vehicle,  among  the  most  significant  being  (Refs.  3,  U): 

(1)  gravitational  torques  due  to  gravity  gradients  across  the  vehicle, 

(2)  aerodynamic  torques  due  to  motion  throufii  the  rarefied  atmosphere, 

(3)  electromagnetic  torques’ due  to  motion  throu§^  a  magnetic  field  and 
surface  charge  on  the  vehicle,  (U)  "radiation  pressure"  torques  due  to 
motion  in  a  photon  stream  from  some  body  such  as  the  Sun,  and  (5)  mechan¬ 
ical  torques  due  to  moving  parts  internal  to  the  vehicle. 

While  it  is  generally  true  that  most,  if  not  all,  of  these  torques 
would  be  insignificant  at  the  Earth's  suface,  they  predcaninate  on  an 
orbital  vehicle  owing  to  the  absence  of  all  other  major  influences.  In 
the  following  sections,  a  discussion  of  these  torques  is  presented. 

2.2  Gravitational  Torques 


The  influence  of  the  gravity  gradient  is  perhaps  best  evidenced  by 
the  Earth's  natural  satellite,  the  Moon,  which,  due  to  an  asymmetrical 
mass  distribution,  always  keeps  approximately  the  same  side  pointed  toward 
the  Earth.  The  recent  literature  is  rich  in  material  on  the  use  of  these 
torques  for  satellite  attitude  control  (Refs.  3-10),  although  as  mentioned 
previously,  much  of  the  work  must  be  classed  as  feasibility  studies  since 
overly  restrictive  assumptions  (such  as  omission  of  all  other  disturbing 
torques,  circular  orbit,  etc.)  have  been  made.  Since  the  Earth's  gravita¬ 
tional  field  may  be  utilized  for  these  purposes  even  at  great  heights, 
however,  it  is  not  surprising  that  the  problem  has  received  so  much 
attention. 

The  general  development  of  the  gravity  gradient  torques  which  act 
on  an  orbiting  vehicle  has  been  performed  by  Roberson  (Refs.  9,  10) 
considering  the  gravitational  potential  which  corresponds  to  the  oblate 
Earth,  but  his  discussion  has  been  unnecessarily  complicated  by  the  intro¬ 
duction  of  numerous  Euler  angles  which  are  particular  to  the  problem  of 
satellite  vehicle  orientation.  DeBra  (Ref.  12),  in  essentially  repeating 
this  development  for  a  spherical  Earth,  arrives  at  an  incorrect  res\ilt 
irtxich,  in  the  limit  of  small  deviation  angles,  merges  with  the  correct 
expressions  given  by  Roberson.  More  recently,  Nidey  (Ref.  20)  has 
published  a  concise  discussion  of  the  gravitational  torques  acting  on  a 
rigid  body. 

For  coD;}leteness,  the  e:q)ressions  for  the  gravitational  torque 
components  along  arbitrary  axes  fixed  in  a  rigid  body  are  derived  from 
first  principles  in  Appendix  A.  Although  only  a  simple  inverse -square 
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attracting  force  field  is  considered,  the  extension  to  include  oblateness 
effects  or  other  anomalies  in  the  force  field  may  be  readily  performed. 
The  result,  as  used  in  Section  1.3,  Eq.  (12),  is 


L»a 

<r^ 

_J 

1 

^Ch-'h) 

vhere  the  direction  cosines  between  the  radius  vector  r  and  the  principal 
body  axes,  r»/v-,  r.j,/r  and  are  given  by  Eqs.  (13),  It  m^  be 

noted  that  when  &  ,  t/-  ,  ^  are  small  amgles  such  that  sin  V  “ 
and  cos  V  ^  I  and  squares  and  products  of  the  small  angles  are  negligible, 
Eqs.  (12)  and  (13)  combine  to  yield. 


L- 

rs  J.  iilL 

(l  -  I  )  Sirt 

z  r* 

>> 

o 

Thus,  as  indicated  previously,  there  is  no  gravitational  yawing  torque  for 
small  angular  deviations  of  the  body  from  the  orbital  coordinates. 

To  illustrate  the  behavior  of  the  gravitational  torques,  ass\xme  that 
^  S  o  and  consider  the  torque  about  the  ^  -axis,  which  is  normal  to  the 
orbit  plane.  The  appropriate  geometry  is  shown  in  Fig.  5.  In  this  case, 
Eq.  (2li)  becomes 


which  is  shown  sketched  qualitatively  in  Fig.  $. 

Because  depends  upon  sinE9  ,  the  gravitational  torque 

vanishes  at  ©»  o  ,  ir/t  ,  tt  ,  •••  .  Accordingly,  these  represent 
possible  equilibrium  positions  of  the  vehicle.  Note  further,  however. 
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that  only  for  <o  is  the  equilibrium  stable,  since  in  this 

case  a  small  positive  increase  in  ©  yields  a  negative  (or  restoring) 
torque  on  the  body .  Thus ,  when  ,  ©»o  ,_Tr  z-TT  ,  •  •  • 

represent  positions  of  stable  equilibrium  while  for 

OT/i  ,  •••  represent  stable  positions.  The  gravitational  torques,  there¬ 
fore,  imply  a  position  of  stable  equilibrium  only  when  the  body's  axis  of 
least  inertia  points  tovrard  the  force  center  (Refs.  3,  h) . 


Figure  5.  Behavior  of  Gravity  Gradient  Torque 


The  motion  of  a  rigid,  dumbbell-shaped  satellite  under  the  influence  of 
the  gravity  gradient  torques  alone  has  been  studied  extensively  in  the 
literature;  for  example,  Schindler  (Ref.  5),  Stocker  amd  Vachino  (Ref.  7) 
and  Baker  (Ref.  8).  It  is  of  interest,  therefore,  to  note  how  the 
present  analysis  may  be  reconciled  with  these  previous  investigations. 

Referring  to  Fig.  6  for  the  geometry  of  the  dumbbell-shaped  satel¬ 
lite,  it  is  clear  that 

Ij  =  7  -  (C 
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hence  for  l>  ^ 
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2  ^ 
2  rs 


Sm  2© 
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Stn  2© 


Comparing  this  result  with  the  second  of  Eqs.  (22),  it  is  seen  that  in 
the  subsequent  work,  the  pitching  motion  of  a  dumbbell-shaped  satellite 
may  be  recovered  by  setting  r^\=.  i  and  neglecting  all  aerodynamic 
influences. 


Figure  6.  Geometry  of  Dumbbell -Shaped  Satellite 


2.3  Aerodynamic  Torques 
2.3.1  General  Considerations 


The  open  literature  contains  ^rprisingly  little  detailed  work  on 
the  influence  of  the  aerodynamic  torques  on  the  vehicle's  motion,  in 
spite  of  the  fact  that  for  Earth  these  torques  are  known  to  be  significant 
over  a  wide  range  of  altitudes  (c.f.  Refs.  U-lh).  Because  of  highly 
rarefied  nature  of  the  atmosphere  at  altitudes  where  a  satellite  may  make 
one  or  more  successful  orbits,  the  molecular  mean  free  path  is  generally 
much  larger  than  any  characteristic  body  dimension  and  free-molecular 
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flow  exists.  Aero(Jynamic  forces  and  moments  on  vehicles  in  this  flow 
regime  are  due  solely  to  the  molecular  bombardment  and  subsequent  reflec- 
ticffi  from  the  siorface  and  are  most  conveniently  predicted  within  the 
framework  of  the  kinetic  theory  of  gases.  The  'fexact" theory  has  been 
developed  extensively  by  numerous  authors  but  the  expressions  for  the 
free-molecular  pressure  and  shearing  stress  coefficients  on  an  arbitrary 
element  of  body  surface  are  generally  so  complicated  that  the  integration 
to  obtain  forces  and  moments*  has  been  performed  for  only  the  simplest 
bodies  (c.f.  Ref.  1,  Section  1). 

More  recently,  an  approximate  free  molecular  theory  has  been 
published  (Ref.  21)  which  permits  calculation  of  the  conventional  aero- 
cfynamic  stability  derivatives  of  arbitrary  bodies  of  revolution.  Davison, 
in  Volume  II  of  this  series  of  report.s  (Ref.  1)  presents  a  detailed 
application  of  this  and  similar  methods  to  a  wide  variety  of  body  shapes. 
The  results  of  these  studies  are  summarized  briefly  in  Section  3.3  of 
this  report. 

2.3.2  Orientation  of  the  Relative  Wind  Vector 

The  aerodynamic  torques  acting  on  a  satellite  may  be  determined 
onge  the  relative  wind  vector  has  been  oriented  with  respect  to  the 
vehicles.  The  relative  wind  vector  is  the  velocity  of  the  vehicle  with 
respect  to  the  local  atmosphere,  and  due  to  atmospheric  motion  relative  to 
inertial  space,  this  vector  may  differ  appreciably  from  the  inertial 
velocity  vector  in  direction  as  well  as  in  magnitude.  Although  a  number 
of  possible  alternatives  suggest  themselves,  the  conventional  angle  of 
attack,  f  and  the  sideslip  angle,  P  ,  are  most  convenient  for  orient¬ 
ing  the  relative  wind  with  respect  to  the  bocfy  axes.  By  retaining  these 
angles,  the  aerodynamic  moment  coefficients  may  be  carried  over  directly 
from  conventional  aerodynamic  theory  and  the  interpretation  of  the 
various  terms  idiich  arise  is  considerably  simplified. 

Referring  to  Fig.  7,  the  components  of  the  relative  velocity 
vector  Vn  may  be  expressed  in  the  body  axes  by  the  matrix  equation 


’v 

Co-0  Co>e  ^ 

Vr 

Stn^ 

SiKt  o6  C»«  ^ 

(25) 


Figure  7.  Orientation  of  the  Relative  Wind  Vector 
with  Respect  to  the  Vehicle 


The  relative  velocity  components  in  the  body  axes  may  also  be 
obtained  through  an  independent  transformation  which  involves  the  previ¬ 
ously  introduced  orientation  angles  &  ,  .  To  obtain  this 

relationship,  first  orient  the  relative  wind  vector  with  respect  to  the 
orbital  coordinates  by  two  new  angles  T  and  A  ,  as  shown  in  Fig.  8. 
These  angles  are  completely  specified  by  the  motion  of  the  vehicle's 
mass  center. 

In  this  case  the  components  of  in  the  orbital  coordinates 

may  be  given  in  matrix  form  as 
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Figure  8.  Orientation  of  the  Relative  Wind  Vector  with  Respect 
to  the  Orbital  Coordinates 


Now,  introducing  the  transformation  from  the  orbital  coordinates  to  the 
body-fixed  axes  as  given  by  Eq.  (b),  there  results 


ci-C(e-'l^)CA  +  S-fsA 

+  c^sVc(©-'!r)cA  ~c4ci'SA 

C^S^©-y)CA-  Sj^S-vAcC6>-y)CA->-  S(^CV-SA. 

Equating  corresponding  components  of  Eqs.  (2?)  and  (27)  gives  the  desired 
relationships  for  Ql,  and  ^  in  terms  of  other  known  angles.  Unfortun¬ 
ately,  the  complexity  of  these  relationships  for  the  general  large  angle 
case  virtually  precludes  their  use  in  analytic  work.  Accordingly,  it  is 
desirable  to  introduce  the  assumption  of  small  angles  consistent  with  idiat 
was  done  in  Section  l.l^.  In  order  to  do  this,  it  must  first  be  shown 
that  T  and  A  are  truly  small. 

The  atmosphere  of  Earth  essentially  rotates  as  a  rigid  body  attached 
to  the  Earth.  Although  it  is  probably  true  that  at  sufficiently  great 
distances  such  motion  will  not  occur,  at  altitudes  >diere  the  aerodynamic 
influences  predominate  the  assumption  of  rigid  body  atmospheric  rotation 
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should  be  acceptable.  The  presence  of  local  vrind  anomalies  may  be 
neglected,  since  they  are  generally  fluctuating  in  magnitude  and  direction 
and  are  not  sufficiently  well  defined  at  satellite  altitudes  to  allow  a 
satisfactory  generalization  to  be  made.  Accordingly,  the  relative  velocity 
vector  may  be  related  to  the  inertial  velocity  vector  through  the  equation 

=  V  -  W  lOg  xT  (28) 

where  cOg  is  the  angular  velocity  vector  of  the  Earth,  V  is' the 
inertial  velocity  vector  and  k  is  a  coefficient  which  may  depend  upon 
altitude  and  permits  deviation  of  the  atmospheric  motion  from  a  purely 
rigid  body  rotation.  The  components  of  ZOg  in  the  orbital  coordinates 
are  found  from  Eq.  (2)  to  be 


=  -COg 

1 

(29) 


while  by  definition  of  the  orbital  coordinates,  r  lies  along  -  55,  . 

The  vector  product  uTgx'r  given  in  Eq.  (2)  then  becomes,  in  matrix  form, 


[^cog  X  =  x~ 
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(30) 


•  •  • 

Now  previous  arguments  have  shown  c5l  ,  i  and  «o  to  be  negligioly 
small  and  hence  the  motion  of  the  vehicle’s  mass  center  is  planar.  By 
definition  of  the  orbital  coordinates,  this  plane  is  coincident  with  tne 
Xo-i.  plane,  and  so  v  has  no  component  along  'fo  •  Let  the  components 
of  ^  in  the  orbital  coordinates  be  given  in  matrix  form  by 
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where  Y*  is  the  inclination  angle  of  the  inertial  velocity  vector  (or 
the  tangent  to  the  orbit)  with  respect  to  the  ><,  axis.  Clearly  Y*  is 
also  corapletely  specified  by  the  motion  of  the  mass  center.  From  Krause 
(Ref.  18)  or  Eq.  (B-5  )»  Appendix  B,  there  is  found 


To^  r* 


r"Cr 


e  Stn  V 

I  -t-  e  co-o-vT 


where  e  is  the  orbital  eccentricity.  For  elliptical  orbits  of  small 
eccentricity  (say  less  than  .06),  it  follows  that  Tf *  is  a  small  angle. 

Combining  Eqs.  (30)  and  (31)  with  Eq.  (28),  the  components  of  the 
relative  velocity  vector  in  the  orbital  coordinates  are  found  to  be^ 


—  *“ 

rt03  Co^4  -N/Co-oY^ 

=:  — 

N/  Sinlf* 

and  the  magnitude  of  Vg  is  given  in  terms  of  known  quantities  as 

(^)=  1  -  ^  4  .  (33) 

Sterne  (Ref.  19)  has  also  presented  this  result.  Comparing  Eqs.  (26)  and 
(32)  allows  the  angles  T  and  A  to  be  evgiluated; 


Svn  Y  = 


Sin'V*f  V  \ 

C^A  * 


(3li) 


For  low  altitude  orbits  about  the  Earth,  k  is  near  unity  and  so 
<X>^&  (2  X  103)  while  at  altitudes  for  which  aerodynamic  effects 

are  significant  V  ft’  O  (2  x  10^)  hence  r«Oa/v  is  small  compared  to 
unity.  In  this  case,  the  ratio  V/V^  may  be  developed  from  Eq.  (33)  by 
a  binomial  expansion  in  the  form. 


^It  is  convenient  to  introduce  the  abbreviation  Vcc-Og  denote 

the  angular  velocity  of  atmospheric  rotation. 


4  =  1  4.  +  <7  (35) 

where  the  terms  omitte(i  produce  an  error  of  less  than  1%,  Note  that  the 
last  term  in  Eq.  (33)  contains  the  only  latitude -dependent  contribution 
to  v/v*.!,  and  that  this  term  is  of  order  (rtOa/v)’^  and  hence 
negligible  in  Eq.  (35).  Inserting  this  relation  into  the  first  of 
Eqs  (31;)  1  and  again  neglecting  terms  of  order  there  results 

Sie  A.—  ^  6>i«  4  Co-0  U.  ,  (36) 

This  expression  shows  that  A  is  a  small  angle  so  that  cos  A  ^  1 
Hence  the  second  of  Eqs.  (3U)  becomes,  to  first  order  terms  in  r-iOa/v  , 

SinT  =  I  -t-  co-o'iT*)  , 

Since  it  has  already  been  observed  that  "tP  is  small,  then  it  follows 
that  IP  is  also  small  and  is  given  approximately  by 

l  +  co^l)  .  (37) 

Now  assuming  (/>  ,  3^  ,  ^  ,  'ZT’  and  A  are  all  small  angles,  and 

comparing  Eqs.  (25)  and  (27), 
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^Spherical  trigonometry  yields  the  relation; 
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from  which  it  follows  that  Oi,  and  ^  are  also  small  omd  given  by  the 
linear  equations 


0i,=  ~  TP  ) 

fS,  Si  0^- A.)  . 


(38) 


It  is  important  to  note  that  this  representation  fails  completely  in  the 
case  of  large  angular  deviations. 

2.3.3  Expansion  of  the  Aerodynamic  Torques 

With  the  orientation  of  the  relative  wind  vector  completely 
specified  for  small  angles  by  Eqs.  (38),  the  aerodynamic  moments  may  be 
determined.  Following  the  conventional  practice 
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where  P  is  the  local  atmospheric  density,  *Sir  is  a  reference  area  on 
the  vehicle,  C  is  a  reference  length  on  the  vehicle  and  ,  Cw»  and 

Cy,  are  respectively  the  dimensionless  rolling,  pitching  and  yawing 
moment  coefficients.  It  is  conventional  to  expand  Ojt  ,  .  Cvn  and  Cvi 
about  the  equilibrium  position  as  a  Taylor  series  in  the  variables  OU  , 

^  f  ^  j  T  ,  ^  and  /[  .  When  the  angular  deviation  of  the  bocfy 

from  equilibrium  is  small,  only  the  leading  terms  in  the  series  need  be 
retained,  and  if  the  bo^y  possesses  symmetry,  a  number  of  the  remaining 
coefficients  may  be  shown  to  be  small  or  else  identically  zero.  For  a 
body  of  revolution,  this  "linearized"  representation  takes  the  form, 
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where  the  second  subscript  on  the  coefficient  of  any  terra  indicates 
differentiation  and  evaluation  as  indicated  below: 


I’M 


^  — ►o 


T— o 

Note  that  the  rate  derivatives  have  been  non-dimensionalized  oy  replacing 
the  rate  (i  ,1“  j  etc.)  by  c  /2Vn  times  the  rate  (  6o  c/z.'Vfe.  j 

^<S:/-aV«  ,  etc.)  These  equations  show  a  strong  coupling  between  Cj 

and  Cv\  ,  while  is  essentially  imcoupled  from  them  except  for  the 
term  ^  which  is  generally  small  for  symmetric  configurations.  The 

terms  involving  the  angular  rates  "f  ,  ^  and  fl  represent  aerodynamic 
damping,  which  is  the  only  "natural"  damping  present  in  the  system.  For 
the  purpose  of  simplification,  only  the  following  dominant  terms  shall  be 


retained: 
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since  only  bodies  of  revolution 
are  to  be  considered.  The  negative  si^  must  be  imposed  on  since 

a  stable  body  at  positive  ot  and  ^  (see  Fig.  7)  experiences  a  negative 
pitching  moment  but  a  positive  yawing  moment  and  hence  for  bodies  of 
revolution,  Cyvi*  and  are  equal  in  magnitude  but  differ  in  sign. 

Inserting  Eqs.  (38)  and  (UO)  into  Eq.  (39)  and  expanding  “f  ,  /  and  , 
the  complete  linearized  expression  for  the  aerocfynamic  torques  becomes, 
for  small  angles  and  angular  rates. 
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where  sin  A  —  A  is  given  by  Eq.  (36),  'jT  is  given  by  Eq.  (37)  and  the 
magnitude  of  is  given  to  a  first  order  in  r-uo^/s/  from  Eq.  (35)  by 

sv  ^  ^  (1,2) 

In  the  following  work  the  coefficients  Of  ,  and  g  , 

which  are  generally  functions  of  altitude  and  velocity,  will  be  regarded 
as  constants  of  the  motion.  The  error  introduced  in  this  way  is  small 
compared  to  the  uncertainties  in  other  quantities  such  as  the  atmospheric 
density. 

2.h  Other  Disturbing  Torques 

2.U.1  Solar  Radiation  Pressure  Torques 

Following  the  gravity  gradient  and  aerodynamic  torques,  the  next 
most  significant  external  torque  appears  to  be  that  due  to  the  solar 
radiation  pressure.  For  extremely  high  altitude  orbits,  this  torque  will 
doubtless  transcend  all  other  predictable  external  torques.  In  this 
connection,  the  work  of  Frye  and  Stearns  (Ref.  U),  Sohn  (Ref.  15),  and, 
more  recently,  Newton  (Ref.  16),  provide  suitable  documentation.  Accord¬ 
ingly,  the  significant  question  in  regard  to  the  present  study  concerns 
the  altitudes  at  which  the  solar  radiation  pressure  will  induce  torques  of 
comparable  magnitude  to  the  aerodynamic  and/or  gravitational  torques. 

Since  the  solar  torques  arise  from  the  momentum  flux  imparted  to 
the  vehicle  by  the  photon  stream  emanating  from  the  Sun,  they  are  amenable 
to  analysis  using  suitably  modified  free  molecular  aerodynamic  concepts. 

A  detailed  examination  of  these  torques  for  several  classes  of  elementary 
boc^  shapes  is  presented  in  Appendix  A  of  Ref.  1.  Fig.  9  summarizes  the 
results  of  these  calculations  by  comparing  the  magnitudes  of  the  initial 
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torque  slopes  liue  to  aerodjmamlc  and  solar  influences  for  three  optimized 
body  shapes^.  It  is  seen  that  for  these  shapes «  the  solar  torques  do  not 
become  con;>arable  to  those  due  to  aerodynamic  effects  until  altitudes  in 
excess  of  UOO  miles  are  achieved.  Wall  (Ref.  11)  in  comparing  directly 
the  magnitudes  of  the  solar  radiation  pressure  and  the  dynamic  pressure 
due  to  the  satellite's  motion  through  the  atmosphere^  also  concludes  that 
aerodynamically  stable  bodies  will  "trim”  into  the  Sun  at  altitudes  above 
about  iiOO  miles. 

Because  the  solar  photon  stream  is  not  conpletely  reflected  by  the 
body>  one  component  of  the  solar  flux  appears  as  a  shearing  stress  on  the 
bocfy-  surface.  This,  in  turn,  results  in  a  residual  torque  arising  from 
surface  tractions  alone  which  threaten  to  perturb  all  but  spherically 
shaped  satellites  vdth  mass  center  at  the  geometrical  center  of  the  body. 
Consequently,  the  solar  torques  may  be  expected  to  be  a  significant  factor 
even  with  a  gravity  gradient  stabilized  vehicle.  This  point  is  illustrated 
in  Fig.  10  where  the  maximum  values  of  the  torques  due  to  aerodynamic, 
gravitational  and  solar  influences  are  compared  for  the  copical  body  shape. 
It  is  seen  that  the  solar  and  gravitational  torques  are  of  conparable 
magnitude  and  significantly  smaller  than  the  aerodynamic  torques  for  alti¬ 
tudes  below  about  IjOO  miles.  Above  600  miles  altitude,  the  atmospheric 
density  becomes  so  small  that  the  gravitational  and  solar  torques  dominate 
the  vehicle's  motioi.  Although  the  absolute  torque  levels  change  with  the 
geometry  of  the  body,  these  qualitative  relationships  remain  true  for  other 
configurations  (Ref.  1).  It  appears,  therefore,  that  the  combined  influ¬ 
ence  of  solar  and  gravitational  torques  (rather  than  either  separately) 
are  required  to  properly  analyze  the  passive  angular  response  of  a  high- 
altitude  satellite. 

As  a  result  of  the  foregoing,  subsequent  work  in  this  study  is 
limited  to  altitudes  below  UOO  miles.  The  ftamework  which  is  developed, 
however,  is  readily  adaptable  to  the  analysis  of  high  altitude  orbits. 


^For  these  calculations,  the  incident  solar  radiation  pressure  is  taken 
as  10“7  psf  and  the  reflectivity  of  the  body  surface  material  is  assumed 
to  be  0.5. 


Absolute  Value  of  Initial  Torque  lope 


Altitude,  h  ,  ^  Statute  Milee 

Figure  9.  Comparison  of  Solar  and  Aero<trnamlo  Initial- 
Torque  Slopes 
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Figure  10.  Maximum  Values  of  Various  Torques 
on  a  Conical  Satellite 


2.U.2  Additional  Torques 


Compared  with  the  aeroc^namiC)  gravitational  and  solar  torques,  the 
remaining  torques  of  a  "predictable"  nature  are  generally  conceded  to  be 
very  small  (Ref.  3U).  For  a  rapidly  spinning  vehicle,  the  most  significant 
additional  torque  appears  to  be  due  to  the  interaction  of  eddy  currents 
induced  in  the  conducting  skin  of  the  satellite  with  the  Earth's  magnetic 
field.  The  problem  has  not  been  treated  extensively  (e.g.,  Vinti,  Refs.  22, 
23,  and  Roberson,  Ref.  3)  to  date  and  appeaurs  to  require  further  stuc(y. 

An  estimate  of  the  spin  retardation  produced  by  these  electromagnetic 
effects  on  Sputnik  III  (Ref.  31)  indicates  this  torque  to  be  roughly 
comparable  to  the  aerodynamic  damping-in-roll  and,  as  will  be  seen  in  the 
next  section,  this  is  an  exceedingly  small  quantity.  For  the  non-spinning 
vehicles  below  hOO  miles  altitude  vAiich  are  considered  here,  electro¬ 
magnetic  effects  appear  insignificant. 

There  are  numerous  other  torques  of  a  statistical  or  random  nature 
(e.g.,  micrometeorites  impacts,  solar  flares,  internal  moving  parts,  etc.) 
which  cannot  be  readily  predicted  and  which  will,  therefore,  be  neglected 
a  priori  in  this  study.  These  omissions  should  not  significantly  affect 
the  subsequent  results.  It  should  be  noted,  however,  that  the  use  of 
internal  moving  parts  may  well  be  among  the  most  efficient  means  of 
achieving  long  duration,  precise  attitude  control.  In  this  connection, 
the  rather  thorough  work  of  Roberson  (Ref.  2h)  is  noted.  The  subject  of 
active  control  system  synthesis  has  not  been  specifically  considered  in 
this  study. 
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3.  EQUATIONS  OF  MOTION  FOR  SLIGHTLY  ECCENTRIC  ORBITS 


3.1  Small  Angle  Equations  for  Slightly  Eccentric  Orbits 

Since  only  the  aerodynamic  and  gravitational  torques  are  to  be 
considered  in  the  subsequent  analysis,  Eq.  (Ul)  may  be  substituted  into 
Eqs.  (22),  for  the  case  of  a  non-spinning,  axisymmetric  vehicle  with  small 
angular  deviations,  to  obtain. 


’(j>  =  pve  r  i 

4  1*  J 


e  .  5^  (e 

L  (U3) 

+  u  -t- 1 )  9^ 

which  are  the  equations  to  be  solved.  Again  it  is  to  be  noted  that 
these  equations  are  linear.  As  noted  previously  for  small  angles  the 
equation  for  the  pitching  motion  is  uncoupled  from  the  equations  for  the 
rolling  and  yawing  motions  and  may  be  solved  sepaurately. 

Now  if  the  secular  effects  of  aerodynamic  drag  on  the  argument  of 
perigee  tO  are  sufficiently  small  (see  Section  1.2)  so  that  to  S  o  , 
then  u.  —  v*  and  ii  =  v  .  It  is,  therefore,  convenient  to  transform 
the  independeiit  variable  from  time  to  the  true  anomaly,  XT  .  Since  the 
motion  of  the  vehicle's  mass  center  is  assumed  to  be  specified  and  known, 
'tr  is  a  specified  function  of  time  alone.  The  necessary  transformation 
relations  are 
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dv 


dt 


In  addition,  it  is  convenient  to  introduce  the  abbreviations 
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whereupon,  applying  Eqs.  (UU)  and  (lj5)  to  Eqs.  (U3)!  there  results 
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where  primes  denote  differentiation  with  respect  to  V' ,  In  order  to 
proceed,  it  is  necessary  to  express  the  coefficients  of  Eqs.  (U6)  as 
functions  of  ir  and  other  constants  of  the  motion.  Owing  to  the  assump¬ 
tion  that  the  center  of  mass  motion  is  uncoupled  from  the  angular  motions 
about  the  mass  center,  the  usual  Keplerian  orbit  elements  may  be  employed. 
The  details  are  tedious  but  straightforward  and  may  be  foiind  in  Appendix  B. 
The  results  are 
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where  the  subscript  P  denotes  that  the  quantity  is  evaluated  at  perigee, 
€  is  the  orbital  eccentricity  and  CP^  ,  g  and  are  constants 
defined  by,  ^ 

St  C 


CP-  s  -  p.  r-^ 

p  ’*’ot  2  1  ^ 

2CR,® -% . 
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The  ratio  of  relative  to  inertial  velocity  Vn,/>/  may  be  evaluated  from 
Eq.  (Ii2)  and  Eqs.  (B-8)  and  (B-9)  of  Appendix  B  as 


1- (W) 

V  -VTj,  (_l  4-fe  C-.roV'J  , 

The  above  expressions  include  only  the  assumption  that  tOa/vp  is  small 
compared  to  unity  and  that  V  ,  V*  and  A.  are  small  angles.  Unfortun¬ 
ately,  Eqs.  (Ii7)  are  much  too  conplicated  for  analytical  work,  and  must  be 
further  simplified.  The  appropriate  simplifying  assumption  is  that  the 
orbital  eccentricity,  6  ,  is  sufficiently  small  (say  e  4  o.i  )  so  that 
e"*-  «■  I  .  This  assumption  is  consistent  with  the  small  angle  restriction 
made  earlier.  Neglecting  squares  of  6  and  j^oducts  of  the  form  ea?a/v- 
the  above  coefficients  simplify  to  the  following  first  order  form: 
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and 
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Note  that  to  this  order  of  approximation,  Vr/v  is  constant  and  equal  to 
its  value  at  perigee.  The  variation  of  d  as  a  function  of  perigee  alti¬ 
tude  and  orbit  inclination  is  shovm  in  Fig.  11. 

In  subsequent  discussions  relating  to  the  numeric  integration  of  the 
"exact”  equation  of  the  pitching  motion,  it  is  the  second  of  Eqs.  (U6) 
with  coefficients  given  by  Eqs.  (50)  and  (51)  vdiich  is  solved.  It  has 
been  shown  in  Section  1  that  this  equation  exactly  describes  the  small 
angle  pitching  motion  of  a  non-splrjiing,  axisymmetric  vehicle  when  the 
rolling  and  yawing  motions  are  suppressed.  Moreover,  the  restriction  to 
small  orbital  eccentricities  poses  no  serious  limitation  on  the  results. 
This  equation  will,  therefore,  be  regarded  as  "exact"  in  relationship  to 
the  approximate  results  to  be  obtained  in  the  following  sections. 

3.2  Approximate  Atmospheric  Density  Variation 

In  order  to  complete  the  evaluation  of  these  coefficients,  it  is 
necessary  to  evaluate  f*  /  (p  .  Since  the  Earth  is  an  oblate  spheroid,  the 
Earth's  atmosphere  is  also  oblate  and  hence  a  constant  density  surface  in 
the  atmosphere  is  not  a  sphere.  Although  this  latitude  variation  of 
density  is  partially  accounted  for  by  representing  =  fC'T-'Te) 
where  rg  varies  with  latitude,  two  other  equally  significant  variations, 
imposed  by  the  latitude  variation  of  surface  gravity  and  temperature,  are 
omitted.  In  order  to  retain  the  maximum  simplicity,  it  will  be  assumed 
here  that  the  Earth  is  spherical,  and  hence  f(U)  where  V\=.  r--a£. 

It  is  now  well-known  that  the  atmospheric  density  variation  above 
about  UOOjOOO  foot  altitude  is  best  approximated  by  a  power  function  of 
the  altitude  (Ref.  25).  However,  under  the  assumption  of  small  orbital 
eccentricity,  it  is  also  possible  to  approximate  the  density  variation 
from  perigee  to  apogee  by  an  exponential  relation  of  the  form 


e 


(52) 


where  C!a  and  are  constants  fixed  by  the  perigee  altitude.  The 
critical  condition  for  Eq.  (52)  to  be  a  satisfactory  representation  is 
that  its  error  be  no  more  than  the  uncertainty  in  the  density  itself 
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over  the  altitude  range  of  interest.  Based  upon  Minzner's  1959  model 
atmosphere  (Ref.  26)  the  uncertainty  due  to  seasonal  and  daily  variations, 
as  well  as  uncertainties  due  to  position  on  the  Earth  is  at  least  a  factor 
of  2  at  these  altitudes.  Above  about  100  statute  miles  altitude  and  for 
orbital  eccentricities  less  than  about  .05,  the  constants  of  Eq.  (52)  m^y 
be  found  such  that  the  error  •’s  within  these  limits  and  hence  the  use  of 
this  representation  is  justified.  The  appropriate  values  of  are 

shown  in  Fig.  12. 

From  the  equation  of  the  elliptical  orbit  (Eq.  (B-2)  of  Appendix  B 
or  Ref.  18) 
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\  +  €  Coo-X/- 


hence 


which,  to  terms  of  order  6  is 
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(52),  and  defining  \  ^  y 
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Now  note  that  when  V*  =  di'n'/x,  c.o-o-tr-=,o  and  the  vehicle  is  at  the  semi- 
latus  rectum,  -p  ,  of  the  elliptical  orbit.  The  altitude  to  this  point  on 
the  orbit  will  be  denoted  here  by  s  T-  3  E  •  Then  from  the  above 
equation. 


Substituting  Eq.  (53)  into  Eq. 
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The  behavior  of  the  dimensionless  atmosphere  param«5ter  \  as  a  function 
of  perigee  altitude  and  orbital  eccentricity  is  shown  in  Fig.  13. 

Substituting  Eqs.  (50)  and  (Sh)  into  Eqs.  (ii6),  there  results  the 
following  equations  for  the  vehicle's  rolling,  pitching  and  yawing  motion: 
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where 
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and,  making  use  of  the  fact  that  Vr/v  is  a  constant  to  a  first  order 
6  as  given  by  Eq.  (5l),  the  constants  (Ji* ,  and  are 

defined  by. 
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Eqs.  (55)  represent  the  fundamental  small  angle  equations  of  motion 
for  an  axisymmetric  satellite  vehicle  moving  in  an  elliptical  orbit  of 
small  eccentricity  and  acted  upon  by  aerodynamic  and  gravitational  torques. 
The  analysis  has  shovm  that  for  a  vehicle  of  this  configuration,  the  motion 
is  characterized  by  six  constants:  the  aerodynamic  parameters  OV*  > 
and  ,  the  inertia  parameter  M  ,  the  orbital  eccentricity  €  and  the 

atmosphere  parameter  X.  •  In  the  following  work,  the  nature  of  the 
vehicle's  angular  motion  will  be  studied  in  terms  of  these  significant 
pararaeters  by  seeking  analytical  solutions  to  Eqs.  (55).  In  the  case  of 
the  general  elliptical  orbit,  only  the  equation  for  the  pitching  motion 
will  be  solved,  but  before  entering  into  this  solution,  it  is  instructive 
to  examine  the  complete  motion  under  cerjtain  simplifying  assumptions. 
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Figure  12.  Varletlon  of  Exponential  Factor,  X.^  ,  with  Perigee  Altitude 
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Figure  13.  Ataosidiere  Parameter,  X  ,  as  a  Functicn  of  Perigee 
Altitude  and  Orbital  Eocentrlcity 


3.3  Order  of  Magnitude  of  Aerodynamic  Terms 


Since  Eqs.  (55)  will  be  employed  for  the  subsequent  analytical  work, 
it  is  desirable  to  effect  the  maximum  simplification  possible.  This 
simplification  can  be  performed  by  identifying  those  terms  remaining  in 
Eqs.  (55)  which  are  as  small  or  smaller  than  the  second-order  squares  and 
products  of  small  angles  previously  neglected  in  the  analysis. 

The  order  of  magnitude  of  the  three  aerodynamic  parameters  P  ,  G2. 
and  R  may  be  estimated  from  Eqs.  (56)  as 


P  =  OfCPp) 

^  -  6) 

P.  =  o  ( C)2p  -  e) 

where  oc  )  is  to  be  read  "the  order  of  magnitude  of  the  enclosed 
quantity"  and  where  /v  has  been  taken  as  unity  for  small  eccentricity 
orbits^ . 

It  is  convenient  to  group  the  vehicle  narameters  in  Eqs.  (U8)  to  form 
new  quantities  which  may  be  more  readi''  ermined.  Define: 


cpp » r 


r  ^ 


_  p  PpTr 


SttC^ 


1  I 


SirC^ 


Ctf  zi: 
I 


Davison  (Ref.  1)  has  calculated  the  attainable  values  of  the  para¬ 
meters  p  and  n  for  an  extensive  number  of  satellite  shapes  using  the 
approximate  free  miolecular  aerodynamic  theory  presented  in  Ref.  (21). 

The  results  of  these  calculations  may  be  summarized  as  shown  in  Figs,  lli 
through  17.  It  may  be  seen  that  values  of  f*  and  in  excess  of 


^Note  from  Eqs.  (5U)  that  (Pi/Pr^^  ^  •  This  quantity  has 

been  set  to  unity  here  to  insure  a  conservative  estimate  throu^out. 
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10^  are  generally  possible  only  for  very  small  vehicles.  Although  values 
of  n,  have  not  been  calculated  specifically,  it  is  reasonable  to  assume 
that  this  parameter  will  not  exceed  the  order  of  Po  .  For  perigee  alti¬ 
tudes  near  or  above  100  statute  miles,  Pp  is  at  most  2  x  10”12  slug/ft^ 
while  Tp  is  of  order  2  x  107  feet.  Employing  these  maximum  estimates 
of  the  various  parameters, 


P  =  OC  PxlO*) 

Q  ^  -6)  (60) 

R  =  ^  ( r^Kio  ^  -  (£) , 


Clearly  P  is  a  potentially  large  quantity  and  must  be  retained. 

Q  and  R  ,  however,  are  seen  to  be  of  order  S  ,  since  the  aerotiynamic 
terms  in  these  quantities  are  at  most  of  order  for  all  reasonable 
satellite  shapes.  Inasmuch  as  this  order  of  magnitude  estimate  has  been 
biased  to  yield  the  largest  values  to  be  expected,  it  appears  conservative 
to  neglect  as  second  order  the  parametric  excitation  terms  (a  R+ 
and  +  2€Str»v)  ^  in  Eqs.  (55)*  The  resulting  equations  then 

become 


+  -»-[3MCl-e<^-v-)+p]©=ZQ  +  €Ps.n-vr'  (61) 

y  -b2Q.y- ^ [[P-m]  ■y-  =■-  P<^  Co^(y+i*:>) 


It  should  be  remarked  that  the  assun?)tion  of  negligible  aeroc^ynamic 
damping  has  been  made  a  priori  in  previous  studies  (e.g.,  DeBra,  Ref.  12 
and  Wall,  Ref.  11),  no  doubt  because  adequate  theoretical  means  to 
predict  the  aerodynamic  damping  characteristics  of  arbitrary  bodies  were 
not  available.  It  is  seen  in  Figs.  1$  and  16,  however,  that  although 
this  damping  is  always  small,  it  need  not  be  negligible  and,  consequently, 
its  first  order  effects  should  be  retained  for  completeness  in  the 
equations  of  motion. 


U8 


Now  introduce  the  exponential  transformations 
,  -  -5RCv)J-w- 

(f>  =  <l>  e.  ^ 

_  —  \  dv- 

0  “  ©  e  •* 

-P  -  I  Q.  t-v)  d-v 


The  typical  derivatives  of  these  functions  are  illustrated  by 

e  =(©  -Q©)e 


"  o'*  -  XQ  ©  '  +  -  Q')^  ^  ^ 


--  ^QC-nt)  d  V 


Inserting  Eqs.  (62)  relations  into  Eqs.  (61)  and  rearranging, 


-  [r^+-r']  ^  -qV'  J  e 


^  £  RCv-i-atv')^^- 


=  ^Z.<5,  +  ^ 


\  Gl,C.v> 

►  •' 

5  Gi,C'«r>4'u- 


V'  •b[^P-M-Q^-<a']-y'=r  P(i  C^,-o(.ir+tO;>  €, 

^CplCvi-ROv)]Ai 

+  (m  +  1)  ~  R  ^  0 

These  transformed  equations  are  particularly  convenient  for  further 
order  of  magnitude  analysis.  A  detailed  examination  of  the  coefficient 
of  0  in  the  pitching  equation  will  serve  to  illustrate  the  procedure. 


It9 


To  terms  of  order  es  »  Eqs.  (56)  yield 


Q 


and 


■z.  2\c^v-^  *  XCo-.\r 


^  ~  Sp  C  {_l  -t- ^€(1  - c»-«v)J  +  2  2p  ^  esio-v-  +  O(e^) 


so  that 


-  /  2\c»«V"  _ 

Q  +Q  =  2p  [l +260- c<«v)J  -  6  C<^^xr- 

“  2*^^  E*  "*  +6^1- c^v)Jsinv~  , 

In  order  for  these  terms  to  be  retained,  it  is  necessary  that  they  be  of 
comparable  magnitude  to  P  .  For  order  of  magnitude  purposes,  as  before, 

Q^  +  q'  =  o(ap  -  €  -  SpX) 

while  from  the  first  of  Eqs.  (58), 


p=  O  (  (Pp)  , 


Now  write 


or, 


-Q--<a'=  P{i-  ^ 
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Introducing  Eqs.  (59),  the  following  order  of  magnitude  estimates 
result: 


Using  the  maximum  values  discussed  previously,  it  is  evident  that  the  first 
two  terms  j'bove  are  negligible  in  comparison  to  unity  for  all  practical 
vehicles  and  orbits.  Because  of  the  appearance  of  P  and  Pp  in  the 
denominator  of  the  last  relation,  however,  this  term  need  not  remain  small 
conpared  to  unity  and  hence  must  be  retained.  Thus,  the  second  of  Eqs.  (63) 
may  be  reduced  to 


— "  r-  /V  T  -  r-  1 

O  -h  €c<«sv +Pj  ©  =  Fz-Q -h  ePsiio-ui0  (61:) 


Applying  the  above  procedure  t£  the  fir^  and  third  of  Eqs.  (63)  and 
neglecting  products  of  the  form  and  R  ^  in  comparison  to  and 

in  the  parametric  excitation,  there  results 


P  +  e  0»-o  -yj-  <p  =  -  ^ 


y-  •+■  1^  P  +  €t<-o^^r■  “MJ  Y'  —  Pd  6 


^QCv-^dli 


(65) 
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It  is  now  necessary  to  consider  the  exponential  terms  remaining  in 
the  right  of  Eqs.  (6U)  and  (6$).  A  detailed  treatment  of  the  effect  of 
exp^S  in  damping  the  satellite’s  pitching  motion  is  presented  in 

Appendix  C,  where  it  is  shown  that  the  time  to  damp  to  one-half  the 
initial  pitch  amplitude  is,  at  most,  comparable  to  the  satellioe's  orbital 
lifetime.  Thus,  the  secular  change  in  cxp(-5Qtv>iw)over  the  time  interval 
for  which  <?p  ,  2?^  ,  C5li»  »  6  etc.,  m^  be  regarded  as  constant  is 

insignificant  and  hence  all  exponential  terms  on  the  right  of  Eqs.  (61i) 
and  (65)  may  be  taken  as  unity  without  serious  error.  Clearly,  if 
^  QCvjiv  changes  only  slightly,  then 

^  [qCv)-  RCv^]«Jv  =  +  -CPO  iv- 


changes  even  less,  since  it  is  proportional  to  the  difference  of  two  small 
quantities  of  approximately  equal  orders  of  magnitude.  The  reduced 
equations  of  motion  then  become. 


+•  €r  “XT  ^ 


©"4-£3M  +p]e  =  i-  ePsinv- 


(66) 


-  //  r- 

-+-LP+€<^ir 


Co-o  (^v+to^  + 


The  entire  effect  of  aerodynamic  damping  is  now  confined  to  the  term 
2.Q  appearing  on  the  right  of  the  pitch  equation,  and  to  the  transforma¬ 
tion  e^ations,  Eqs.  (62),  relating  &  ,  V  and  ^  to  ^ 
and  ^  .  The  apparent  inconsistency  of  having  the  forced  ir.tion  decay 
according  to  Eqs.  (62)  is  removed  trtien  it  is  recalled  that  the  exponentials 
on  the  ri^t  of  Eqs.  (6U)  and  (65)  differ  from  unity  by  terms  proportional 
to  €  ,  and  these  in  turn  give  rise  to  terms  of  order  which  are 
negligible  as  excitation.  If  the  foregoing  arguments  are  applied  to 
Eqs.  (62),  there  results 

6\Si*»'V'Av" 

5  f  ^ 
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which,  for  injection  at  perigee,  becomes 


^  ,  ,  e  C'  “  Co-a  v) 

e  =  »  -  6  +  o(e^) 

9  ip  J 


Thus,  the  only  significant  contribution  of  the  aerodynamic  damping 
is  to  slightly  alter  the  equilibrium  pitch  angle  such  that  the  damping 
moment  (due  to  the  inertial  pitch  rate  of  the  satellite)  is  countered  by 
a  small  static  pitching  moment.  Subsequent  numerical  calculations  based 
upon  the  complete  equation  of  motion  (see  Section  3.1)  amply  confirm  the 
negligibility  of  the  aerodynamic  damping  insofar  as  its  effect  on  decay 
of  the  angular  oscillations  is  concerned,  and  this  may  be  taken  as  a 
posteriori  justification  of  the  assumptions  anc  approximation  introduced 
above . 

The  reduced  equations  of  motion,  Eqs.  (66),  may  be  further  simplified 
without  additional  error  if  the  term  is  removed  from  the  right  of  the 
yaw  equation  as  follows:  Integrate  the  first  of’  Eqs.  (66)  once  with 
respect  to  v- ,  neglecting  the  term  ec-o^v*  .  The  result  is 

^  '  S  +  Ic 

where  k  is  a  constant  of  integration  determined  by  the  initial  condi¬ 
tions.  Substitution  of  this  result  into  the  yaw  equation  leads  to  the 
new  differential  system 


/ 


0"  +  []3M  +Cl-3Vl)e<^v- ep&m-u-  (67) 

'Z"  '  +  6  -V  =■  Pci  Co-o  C'Lr+uo')  +kCM-1-l)  . 

These  equations  are  now  uncoupled  sufficiently  to  permit  solution,  since 
©  and  "v^  may  be  determined  from  the  last  two  equat^ns  and  ip' 
then  substituted  into  the  first  in  order  to  solve  for  . 


Note  that  when  the  aerodynamic  influences  are  absent,  the  second 
of  Eqs.  (67)  becomes 

b"  +■  =-zes‘oir- 

vrtiich  is  a  non-homogeneous  form  of  Mathieu’s  differential  equation 
(Refs.  32-35).  That  an  equation  of  this  form  governs  the  small  angle 
pitching  motion  of  a  rigid  satellite  in  an  elliptical  orbit  under  the 
influence  of  the  gravitational  torques  alone  has  been  demonstrated  by 
Baker  (Ref.  8).  This  result  is  also  implicit  in  the  large  angle  equations 
for  the  pitching  motion  derived  Beletskiy  (Ref.  27)  and  by  Frick  and 
Garber  (Ref.  28).  Inasmuch  as  this  equation  depends  only  on  the  two 
parameters  e  and  M  ,  a  unique  solution  may  be  obtained. 

The  independence  of  altitude  of  this  equation  is  due  to  the  transformation 
to  vr  as  the  independent  variable  in  place  of  time. 


SEA  LEVEL  WEIGHT,  W  POUNDS 

Fipjre  llj.  Aerodynamic  Moment  Parameter  T  for  a  Number  of 
Optimum  Bocfy  Shapes,  v/v^  =  >/t.. 


SEA  LEVEL  WEIGHT,  W  --  POUNDS 

Figure  15.  Aerodynamic  Moment  Parameter  P  for  a  Number  of 
Optimum  Body  Shapes,  V/V^  = 
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SEA  LEVEL  WEIGHT,  W  POHMDS 

Figure  16.  Aerodtjmamic  Du^dng  Paruieter  for  a  Number  of 
Optimum  Bodty  Shapes,  V/v^  = '/E. 
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SEA  LEVEL  WKIfflT,  W  POUNDS 

Figure  17.  Aero^jrnamlc  Deaiping  Paraaeter  R  for  a  Nunber  of 
Optlmun  Boc^jr  Shapes,  v/v^s»/4 
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3.U  Fourier  Series  Expansion  of  Periodic  Coefficients 

It  may  be  noted  that  all  of  the  coefficients  of  Eqs.  (67)  are 
periodic  functions  of  the  independent  variable  "vr  vd.th  period  ZIT ,  and 
hence  may  be  expanded  in  Fourier  series.  Moreover,  the  coefficients  of 
the  homogeneous  equations  are  even  functions,  involving  cos  V  only,  and 
thus  the  Fourier  cosine  series  suffices  for  the  expansion  of  these  terms. 
In  order  to  outline  the  procedure,  the  pitching  equation  is  considered 
in  detail. 

Let 


e'‘  ^  ^ 


(68) 


where 


Vcooir 


3M(i -eco-oir)  +€co^v-  +  (ppe  [^i+zeCi- 


and 

=  {ZQ +6  Psinvj 


In  order  to  expand  the  function  Xcv)  ,  assume 

CO 


>diere 


ir 


A.  =  iCCtWt 

J  n 


and 


Evaluating  the  coefficients, 


Ac=  3M  +  (1 +ze)  ^  ^  «it  -26^  ^ 


Xc^t 

S  Co-flt  <it 


and 


An  = 


IT 


(3N1-0\  Co-ot  Co-«v>-t«i't 

•'o 

I 


>  *  \  \  C4M»'t 


Kc«-t 


-  \  e 

IT 


Co-oVlt  dt.  jV)* 


The  integrals  may  be  readily  evaluated  with  the  help  of  the  integral  form 
for  the  modified  Bessel  functions  of  the  first  kind  Bind  order  n  (Ref.  33> 
Page  l5l) 


I„U)  =  \  ^  c^vit  dt  ,  ^  (69) 

The  results  are;^ 

A„  =  3M  +l€[U\)-I,(l)]| 

A,  =  (1  -3m)  6  +2(!J*jl,(l)-€[UX)  -Zl,(b  4  1^0)]] 

A„  =  2CPp»{l„W-€[l„.,(X)-zIl\)  +  I.„a)]]^ 


^The  first  few  modified  Bessel  functions  are  shown  plotted  in  Fig.  18, 
for  values  of  the  argument  up  to  6.0. 


60 


Eqs.  (70)  completely  define  the  in  terms  of  other,  known 

qujintities.  When  €  =  0  ,  X®  ®  ®  circular  orbit  results.  The 

coefficients  in  this  case  all  vanish  with  the  exception  of  A* 

which  has  the  value 


-hC?/  . 


When  X  is  large ,  the  [  may  be  evaluated  'dth  the  help  of  the 
asymptotic  expansion  of  the  modified  Bessel  functions.  From  Wittaker  and 
Watson  (Ref.  3h,  Page  373), 


t 


\»\ 


Inserting  this  re.-ult  in  Eqs.  (70)  and  neglecting  terms  of  order  G/\  , 
there  results 


A,S  3M 

A,  3-  (|-3M)6  +2<SV*I,CX) 

A.  =  2.<Pp^  "• 


One  notes  that  as  €  — ►  o  ,  X~^c>  also  so  that  Eqs.  (71)  have  the 
correct  limiting  values  for  circular  orbits  as  well  as  for  large  X  . 
Thus,  Eqs.  (71)  should  represent  reasonable  approximate  values  for  the 
f over  the  entire  range  of  X  considered.  A  nunerical  comparison 
of  the  exact  and  approximate  coefficients  shows  the  agreement  to  be 
excellent  over  the  entire  range  of  interest,  and  hence  Eqs.  (71)  will  be 
employed  in  the  subsequent  work. 

The  non-homogene ous  portion  of  Eq.  (66),  (v)  ,  must  now  be 

expanded  in  a  Fourier  series; 

yri(y)=  e  -^cpp*  e  _  2|  s«n  ■u-  +  z  . 
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In  this  case,  the  "forcing  function"  is  a  purely  periodic  function  of 
with  period  ZTT  ,  and  may  be  expanded  in  a  Fourier  series  of  the  form, 

eo 

•»?(V)=  {C„S  m  n  V  +  Co-o  n  v"  j 


The  coefficients  are  evaluated  in  the  same  manner  as  before: 

c,=  Z6{^'r,a) 

Do  =  z  ap*i<.(A) 

D„  =  4  ,  11=1,1,5,  •••  _ 

The  complete  pitching  equation  thus  assumes  the  form, 

e  +  |[3M  +  +2CP*  l,()0]co^-vr+ ...Je 

+  Sinv-  +  ..  • 


n  =■  z,5j4j*  •• 

(72) 


The  yawing  equation,  given  by  the  third  of  Eqs.  (6?)  is  amenable  to 
a  similar  ainalysis,  the  result  being 

|^[i  -»-C?p*I^(X)]  +  [iC5?.*I,CX)  V' 


=  k(M+-i)+  20V*^ci  Ilex')  co^io 


(7U) 


+  <3p*d[l,CX)  +1^0)]  c<M»«Oc*<.\r'  -2,^  d  X,U) 
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Figrare  18 ,  Kodified 


Bessel  Functions  Versus  X 
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U .  ANGULAR  MOTIONS  ON  A  CIRCULAR  ORBIT 


h.l  Existence  and  Stability  of  Equilibrium  Positions 


The  general  equations  of  motion  developed  in  the  last  section  apply 
to  motion  along  elliptical  orbits  of  small  eccentricity  and  include  the 
first  order  effects  of  aerodynamic  and  gravitational  torques  and  aero¬ 
dynamic  danping.  These  equations,  however,  possess  "simple”  solutions 
only  in  exceptional  circumstances  and  even  in  these  cases,  considerable 
analysis  is  required  before  the  complete  vehicle  motion  may  be  determined. 

In  order  to  gain  some  insight  into  the  complete  angular  motion  of 
an  ar+ificial  satellite  without  such  elaborate  analysis,  it  is  instructive 
to  seex  solutions  of  Eqs.  (6?)  for  the  case  of  a  circular  orbit 

V.  =  o  .  The  assumption  of  a  circular  orbit  has  a  certain  amount 
of  practical  significance  since  such  orbits  are  desirable  for  numerous 
satellite  applications  (e.g.,  weather  prediction,  military  reconnaissance, 
etc.)  where  stabilization  of  the  vehicle  is  required.  It  must  be  noted, 
however,  that  an  exactly  circular  orbit  is  impossible  considering  the 
perturbations  due  to  planetary  oblateness  and  atmospheric  drag.  As  a 
result,  this  assumption  must  be  justified  primarily  on  the  degree  of 
mathematical  simplification  tdiich  it  affords. 

Under  the  assumptions  of  a  circular  orbit,  Eqs.  (67)  sinqjlify  to 

read 


=  Z  a*  (75) 

=  C5V*<i  4-  kC^+O 

where  the  first  of  these  results  form  a  direct  integration  of  (p  ®  > 

as  before.  Clearly,  if  k  is  non-zero,  then  the  rolling  velocity, 
contains  a  constant  (secular)  terra  which  prevents  ^  from  remaining  small. 
In  this  sense,  the  first  of  Eqs.  (69)  establishes  an  approximate  restric¬ 
tion  upon  the  initial  conditions  which  must  be  satisfied  if  the  i^ll  angle, 
^  ,  is  to  remain  indefinitely  bounded;  nivT.ely,  that  ?  fo)  ■=  -  Tpfo)  . 
Since  practical  accuracy  requirements  will  probably  prevent  this  precise 
condition  from  being  met,  it  appears  that  an  auxiliary  roll-stabilization 
system  will  be  necessary  if  it  is  required  that  a  body  of  revolution  orbit 
with  one  side  directed  toward  the  planet. 
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Two  imfjortant  properties  of  the  motion  may  be  deteriid.ned  directly 
from  Eqs.  even  before  a  complete  solution  is  obtained;  these  are 

the  existence  and  stability  of  equilibrium  orientations  of  the  satellite. 
As  mentioned  above,  the  rolling  motion  grows  secularly  if  W oj hence 
any  equilibrium  solutions  must  necessarily  satisfy  ii^ich 

implies  y  *  o  and,  at  most,  ,  a  constant.  From  the  third  of 

Eqs.  (75)  it  is  seen  that  is  a  particular  solution  only  when  (p* 

or  d  vanishes  exactly.  Thus,  equilibrium  solutions  for  the  case  of  an 
aerodynamically  responsive  satellite  exist  only  for  equatorial  orbits  for 
which  d  a  o  . 

The  second  of  Eqs.  (75)  shcn»3  that  an  equilibrixim  pitch  orientation 
is  possible  everywhere  since  S  =  Z  ^®P^®sents  a  pai’ticular 

solution.  The  physical  origin  of  this  term  is  that  a  steady-state  pitch¬ 
ing  moment  must  be  provided  to  counter  the  damping  torque  produced  by  the 
fact  that  even  in  a  circular  orbit,  the  satellite  is  pitching  at  a  rate 
equal  to  its  angular  velocity  in  orbit.  For  practic^  values  of 
and  <P*  ,  however,  this  term  is  insignificant.^ 

Thus,  for  equatorial  orbits,  ccanplete  equilibrium  is  theoretically 
possible.  In  this  case,  the  vehicle  body  axes  remain  essentially  aligned 
with  the  orbital  coordinates  and  the  vehicle  moves  with  the  same  side 
always  directed  toward  the  planet.  Concerning  the  stability  of  this 
equilibrium,  it  is  only  necessary  to  examine  the  coefficients  of  the 
independent  variables  on  the  left-hand  sides  of  Eqs.  (75).  Stability 
exists  if  the  conditions 

Op*  4  >  o 

(76) 

CPp*"  +  »  >  o 

are  satisfied.  It  follows  at  once  that  for  very  slender  shapes  (for 
which  M3?-1  )  pitch  stability  exists  for  C3^*^3  >diile  for  dynamically 
spherical  bodies  ( M  —  ,  pitch  stability  is  indicated  for  all  positive 

values  of  (P*  ,  For  disk-like  bodies,  for  -vdiich  the  inertia  parameter  M 
approaches  unity,  the  axis  becomes  an  axis  of  least  inertia  and  the 
gravitational  torques  alone  are  adequate  to  insure  a  stable  pitch  equili¬ 
brium.  The  yaw  equilibrium  in  circular,  equatorial  orbits  is  stable 
provided  only  that  <Pp*  >  ~l  . 

The  foregoing  stability  considerations  ai'e  shown  graphically  in 
Fig.  19. 


^In  this  connection,  refer  to  Fig.  2h. 


Figure  19.  Stability  Boundaries  of  the  Aerodynamic 
Parameter  for  Circular  Orbits 

p 


The  general  conclusions  may  be  summarized  succinctly  by  noting  that 
for  bodies  of  revolution  with  M  £  73  ,  yaw  stability  is  insured  if  the 
pitching  motion  is  stable,  while  for  m  Z  */3  pitch  stability  is  insured 
if  the  yawing  motion  is  rendered  stable.  In  either  case,  the  rolling 
motion  may  or  may  not  be  bounded  depending  upon  whether  the  precise 
initial  conditions  required  by  the  first  of  Eqs.  (69)  can  be  matched  or  not. 


i 
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U . 2  Approxlnate  Behavior  of  Complete  System 

The  conqslete  analytical  solutions  to  Eqs.  (75)  may  be  readily 
obtained,  since  the  homogeneous  equations  associated  with  the  second  and 
third  of  Eqs.  (75)  are  merely  the  harmonic  equations.  The  complete 
results  are 


0  =  d  co-oA-U"  +•  (B»  S.nA-0-  + 


y-  *  C.0-O  B -vr  4* 

^  ~  SiM  fcv  —  Siw  Cv+wll) 

where 

A^aA«  =  "SM  -V-C5r* 

s  !  -»-(P*  (78) 

k  s  y  ^  -  'i'Co') 

and  ^  ,  <B  ,  <S  ,  &  are  constants  defined  by  the  initial  conditions 
of  the  motion.  Since  the  orbit  is  circular,  no  generality  is  lost  by 
taking  \r=o  as  the  initial  (or  injection)  point,  whereupon  there  results 
(for  A>e>  ,  e.>o  ), 

Q.  -  eco')  +  z  ® 


CB  ^  e'co:>/^ 

C  *=  -  dc»«<0-  •^^C^A+'O  “ 

a  -  isirtw/B  =  -fV©)  /  B 


(79) 


Now  in  general,  injection  into  orbit  will  be  accomplished  in  such  a 
way  that  the  initial  angular  deviations  will  be  small,  as  will  the  rolling 
and  yawing  rates  and  •  The  pitch  rate,  &'co')  ,  however,  will 

generally  be  of  order  unity,  corresponding  to  an  approximate  rectilinear 
injection  path  tangent  to  the  orbit.  These  considerations  suggest  that 
suitable  approximations  to  Eqs.  (77)- (79)  are 
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e  = 


A 


A \r 


V'  ^ 
f-  « 


-  d  s.^  Cv-f-to)  - 


(80) 


The  first  of  Eqs.  (80)  shows  that  although  the  pitching  motion  is 
stable  for  all  such  that  A  io  ,  when  A  becomes  too  small,  the 

pitch  amplitude  may  exceed  that  for  small  oscillations  and  the  entire 
analysis  will  fail.  Likewise,  when  the  natural  frequency  of  the  pitching 
motion,  A  ,  becomes  comparable  to  the  orbital  frequency  (unity)  the 
motion  must  be  considered  effectively  divergent.  Practically  speaking, 
then,  values  of  A  ^  I  will  not  be  of  great  value  in  providing  the  desired 
degree  of  vehicle  orientation,  although  an  estimate  of  the  order  of  these 
effects  will  be  required  to  hi^er  altitudes  even  for  inertially 
stabilized  systems. 


The  terms  on  the  rig^t  of  the  second  of  Eqs.  (7h)  represent  the 
steady-state  motion  of  -vf  .  The  first  term  represents  that  motion  which 
is  induced  by  the  rotatiOTi  of  the  atmosphere.  Over  the  altitude  range  of 
interest,  the  amplitude  of  this  forced  motion  never  exceeds  about  ±li°  as 
shown  in  Fig.  11.  Moreover,  since  this  forced  oscillation  is  periodic  in 
•XT’  with  period  z-tt  ,  it  may  be  removed  by  inertial  devices  such  as 
inertia  wheels,  etc.,  without  the  device  becoming  saturated.  The  second 
terra  represents  a  constant  displacement  of  the  yaw  equilibrium.  Note  that 
when  the  initial  conditions  are  adjusted  so  that  k  •<?  ,  then  the  yawing 
motion  is  essentially  independent  of  the  aerodynamic  characteristics  of 
the  bodyj  that  is,  independent  of  .  The  reason  for  this  independence 

is  that  the  aerodynamic  damping  is  virtually  negligible,  «uid  hence  an 
aerodynandcally  stable  body  aligns  itself  along  the  relative 

wind  vector  and  moves  in  phase  with  it.  When  -i  c<Pp,*<o  ,  the  forced 
yawing  motion  becomes  180°  out  of  phase  with  the  forcing  function. 

Finally,  from  the  third  of  Eqs.  (7li),  it  is  seen  that  the  rolling 
motion  is  90®  out  of  phase  with  the  yawing  motion  and  is  wholly  oscillatory 
provided  that  the  initial  conditions  are  properly  adjusted  so  that  k«-o  . 
If  k  is  different  from  zero,  a  small  secular  term  will  appear,  causing 
^  to  eventually  become  large. 

Thus  when  the  pitching  motion  is  absent,  the  rolling  and  yawing 
motions  may  be  approximately  displayed  on  a  phase  diagram  such  as  shown 
in  Fig.  20(a)  for  k=o  and  as  shown  in  Fig.  20(b)  for  . 
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Figure  20.  Approximate  Behavior  of  the  Forced  Rolling  and 
Yawing  Motions,  (a)  K.=e»  ;  (b)  k^o. 


U.3  Exact  Results  for  Small  Angle  Pitching  Motion 
U.3.1  General 


In  Section  I.U,  it  is  shown  ty  Eq.  (19)  that  if  the  rolling  and 
yawing  motions  are  suppressed,  then  the  pitching  motion  of  a  non-spinning 
axisyinmetric  vehicle  may  be  determined  eaqplicitly  even  for  large  values 
of  ©  .  Since  the  problem  of  obtaining  a  desired  satellite  orientation 
is  principally  concerned  with  the  pitching  motions,  it  is  appropriate  to 
discuss  this  case  in  detail.  Although  only  the  small -angle  equations  are 
discussed  explicitly,  no  fundamental  change  occurs  for  large  angles.  In 
this  connection  it  should  be  noted  that  Beletskiy  (Ref.  27)  has  solved 
this  problem  analytically  in  terms  of  elliptic  integrals  for  the  gravity 
gradient  acting  on  the  dumbbell  satellite  in  a  circular  orbit. 

The  conplete  small-angle  solution  for  the  circular  orbit  pitching 
motion  may  be  found  from  Eqs.  (77)-(79)  as 


1 


(81) 


=  eCcO  C^Air  +  Sin  Av  +X 


where  A*"  and  23*  may  be  expanded  from  Eqs.  (57),  (59)  and  (76)  to 
read 


a'^=3M  ■+  2 

23p*=  Pprp(i-clco*L) 

and  where  P  and  P^  are  defined  by  Eqs.  (59)  as 


(82) 


Stt  ^ 

I 


Thus,  the  vehicle  is  characterized  by  the  three  parameters  T  ,  M  and 
Pj  while  the  orbit  is  specified  by  the  altitude  K  and  the  inclination  C  . 
Since  Eqs.  (8l)  and  (82)  present  the  results  in  such  a  concise  analytical 
form,  the  separate  effects  of  the  various  parameters  may  be  readily  seen. 

For  this  reason,  the  following  discussion  and  related  figures  are  mainly 
for  illustrative  purposes. 


U.3.2  Effects  of  Vehicle  Parameters 


Perhaps  of  greatest  interest  is  the  sunplitude  of  the  pitching  motions 
as  a  function  of  the  vehicle  characteristics  and  the  orbit  altitude. 
Rewriting  Eq.  (81)  in  phase  angle  form,  it  is  readily  found  that 

e,..  =  ®»..=  y 4C®'c»>/a)‘  '  +  2.  VM  . 

Figs.  21  and  22  show  the  results  derived  from  this  equation  for  several 
values  of  M  and  P  ,  with  P©  so  .  Initial  (or  injection)  conditions 
of  SCo'i^o  and  ©''<c>>=:6.  s^are  employed  for  definiteness  but  small 
deviations  from  these  values  will  have  little  influence  on  the  results.  As 
before,  atmospheric  densities  have  been  calculated  from  the  1959  ARDC 
(Minzner)  atmosphere,  (Ref.  26). 
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Fig.  21  shows  clearly  that  even  large  values  of  p  have  relatively 
little  effect  in  containing  the  pitch  oscillations  above  about  UOO  miles 
altitude,  even  for  M>o  .  It  is  seen  that  the  asyiriptotic  (or  P  =•  o  ) 
results  are  very  rapidly  reached  at  or  near  these  heights.  For  those 
shapes  for  idiich  the  gravity  torques  oppose  those  due  to  aerodynamic 
forces  (i.e.,  M  <o  ),  the  deterioration  of  small  angle  motion  occurs  at 
lower  altitudes.  In  all  cases,  however,  it  is  seen  that  for  altitudes 
below  about  l80  miles,  the  effect  of  M  is  negligible,  even  for  the 
relatively  small  values  of  P  .  In  this  sense.  Figs.  21  and  22  permit  the 
definition  of  regions  of  influence  of  the  various  external  torques^: 

Region  I  (below  180-200  miles)  Aerodynamic  torques 
dominate  the  angular  motions. 

\ 

Region  II  (I8O-UOO  miles)  Aerodynamic  and  gravitational 
torques  are  of  comparable  importance. 

Region  III  (1^00-600  miles)  Aerodynamic,  gravitational  and 
solar  torques  are  of  comparable  significance. 

Region  IV  (above  600  miles)  Solar  and  gravitational 
torques  dominate  the  angular  motions. 

Fig.  23  presents  the  typical  effect  of  M  on  the  actual  librations 
of  a  satellite  with  aerodynamic  stability.  As  expected,  the  increased 
stability  associated  with  increasing  M  both  depresses  the  amplitude  and 
shortens  the  period  of  the  oscillations.  For  comparison,  the  corresponding 
curve  for  a  gravity  gradient  stabilized  vehicle  is  also  shown. 

The  effect  of  P j  ,  as  previously  mentioned,  is  too  small  to  be  of 
practical  significance  for  all  reasonable  configurations.  Fig.  2k  presents 
the  equilibrium  angle  z  versus  altitude  for  several  values  of  f’f/P 

in  order  to  substantiate  this  point.  Insofar  as  the  decay  of  the  oscilla¬ 
tions  are  concerned.  Appendix  C  is  to  be  consulted. 

h.3-3  Effects  of  Orbit  Parameters 

The  effects  of  orbit  altitude  are  anqply  displayed  in  Figs.  21  and  22 
and  do  not  require  further  elaboration.  Since  the  orbit  is  circular,  the 
eccentricity  is  absent  and  the  only  remaining  parameter  is  I  ,  the  orbit 
inclination.  Since  i  influences  0  by  slightly  modifying  the  magnitude 
of  the  velocity  of  the  satellite  relative  to  the  atmosphere,  its  effects 


^Specific  control  torques  can,  of  course,  transcend  any  and  all  external 
torques  if  so  desired.  These  have  been  omitted  consistently  in  this  work. 
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may  be  readily  anticipated.  Fig.  2$  presents  a  typical  numeric  result 
for  &  versus  v  to  illustrate  the  effect.  Since  the  order  of 
magnitude  arguments  indicated  secondary  importance  for  the  atmospheric 
rotation,  it  is  not  surprising  that  only  a  small  difference  between  the 
retrograde  (^’  i8o*)  and  progressive  orbits  occur.  In  general, 

it  must  be  concluded  that  for  a  spherically  symmetric  atmosphere,  the 
orbit  inclination  is  of  minor  significance  coirpared  to  (for  exanple)  the 
eccentricity  or  perigee  altitude.  In  the  presence  of  atmospheric  oblate¬ 
ness,  the  "spring  constant"  of  the  system  contains  an  additional  periodic 
contribution  which  considerably  complicates  the  analysis.  Since  the 
governing  equation  in  this  case  is  of  the  Mathieu  type,  the  discussion  of 
Section  5  may  be  used  to  estimate  the  effects  of  atmospheric  oblateness. 


Figure  21.  Combined  Effects  of  Inertia  and  Static  Aeroc^amlc  Moment 
Parameters,  M  and  P  ,  on  Circular  Orbit  Amplitude 
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5.  PITCHING  MOTION  ON  AN  ECCENTRIC  ORBIT 


5.1  Exact  Solution  of  the  Homogeneous  Pitching  Equation 

The  equation  for  the  pitching  motion  on  an  eccentric  orbit,  including 
the  first  order  contributions  from  the  aeroc^amic  and  gravitational 
torques  is  given  Eq.  (73)  as: 


+-<3p^l,CX)]  +[(i  -5M)fe  + 

+  =  2  2^I.(X)  -z6s,^v- 

net-  ^ 

«D 

+  4  2^  ^  S'"  ^  ^  V-]  ^ 

Yl»  I 

The  associated  homogeneous  equation  may  be  brought  into  the  canonical 
form  of  Hill's  differential  equation  (Wittaker  and  Watson,  Ref.  3U, 

Page  l:lii)  through  the  introduction  of  the  transformation 

■f  , 

whereupon  the  homogeneous  equation  becomes 

O  ja,  +  2.  a,  4- z  St  ■+ •  •  •  J  ^  =<p 

where  ao^^Aoi  Sn*  2-An  ,  v»=.|, and  the  primes  now  denote 
differentiation  with  respect  to  c|>  .  With  the  definition  '•'he 

above  equation  may  also  be  written  as 

.  «,  ztncf  -j  _ 

^'’+  {2:  a,e  f©=o  _  (83) 

^  n»<-oo 


This  equation  is  known  to  have  a  solution  of  the  form  (c.f.  Ref.  3h,  Page  i:l6) 


0 


(81:) 
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where  and  all  of  t^ie  are  to  be  determined.  Differentiating 

Eq.  (8ii)  twice  with  respect  to  <p  and  substituting  into  Eq.  (83): 

»  ,  ^+24V;<p  «  •  ^  +  2t(:n^v)]c<> 

e  +££  =0  (85) 

The  operation  of  rearranging  the  series  in  the  second  term  above  is  justi¬ 
fied  since  the  series  given  by  Eq.  (6li)  may  be  shown  to  converge  absolutely^. 
In  order  to  obtain  a  recurrence  relationship  for  the  {  ,  set  "1^=  + 

in  the  second  term  of  Eq.  (85)  and  adjust  the  index  so  that  the  expon¬ 
ential  terms  may  be  factored  out.  Factoring  the  exponential,  setting 
-p  =  V?  and  equating  the  coefficients  of  like  exponentials  yields  the 
recurrence  formula  for  the  {t*,}  ' 


(2 ^  =0^  >;=Oj±.i  5±z 


(86) 


n  =  -eD 


Eqs.  (86)  constitute  an  infinite  set  of  linear  equations  for  the  {t>^} 
and  the  undetermined  exponent  .  Since  the  {3n]  are  completely 
specified,  the  {  may  be  found  onch  is  known.  Wittadcer  and  Watson 
(Ref.  3h)  Page  Ul5)  show  that is  the  root  of  the  equation 

S.n^  =  ^\Co)Sio'^  (4  )  (87) 

where  is  the  infinite  determinant  of  the  coefficients  of 

in  a  set  of  equations  obtained  from  Eqs.  (86)  by  dividing  the  equation 

by  (41^*—  a«)  to  insure  convergence.  is  A{6^)  evaluated  at 

=  o  ,  viz; 


4he  Fourier  series  representing  is  alreaify  known  to  satisfy  this 

condition. 
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(88) 


vhere  So  cannot  be  zero  or  the  square  of  an  even  integer.  is 

called  Hill’s  determinant. 

The  nature  of  the  solution  for  &  ,  given  by  Eq.  (85)  is  grossly 

dependent  upon  the  value  of  — ^  obtained  from  Eq.  (87).  Since  is  a 
root  of  Eq.  (87),  is  also  a  root  and  in  general,  two  linearly 

independent  solutions  to  Eq.  (83)  correspond  to  these  two  roots.  This 
implies  that  if  the  complex  number  contains  a  non-zero  real  part, 

then  one  of  the  two  linearly  independent  solutions  of  the  homogeneous 
equation  contains  a  positive,  real  exponential  factor,  and  hence  this 
solution  is  unbounded  for  increasing  <P  .  Note  that  in  the  cases  where 
the  solutions  corresponding  to  and  are  not  linearly  independent, 

the  second  independent  solution  contains  a  logarithmic  term  and  is  again 
unbounded  (Ref.  32,  Page  556).  Since  the  complete  homogeneous  solution 
contains  a  linear  combination  of  the  two  independent  solutions,  it  is,  in 
principle,  possible  to  find  initial  conditions  for  the  motion  such  that 
the  unbounded  solution  is  absent.  However,  these  cases  are  inherently 
unstable  and  the  complete  solution  is,  therefore,  said  to  be  unstable  for 

56  o  ,  The  problem  thus  becomes  one  of  determining  the  conditions 
under  which  the  solution  of  Eq.  (83)  is  stable;  that  is,  the  conditions  for 
which  is  a  purely  imaginary  nmnber. 

In  order  to  proceed,  define. 
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£50  that,  from  Eq.  (87), 

(89) 

Now  the  {Sw]  are  all  real  so  it  follows  that  is  real  and  hence 

three  cases  are  possible  for  ’§  ,  viz;  (i)  o  ,  (ii)  o<  f  ^<  \ 

(iii)  5'*->  1  .  Eq.  (89)  may  be  readily  solved  for  these  three  cases  to 
yield: 

^  =  ±  S.ntl  i^f)  ;  <0 

lr 

^  ^  ±Zl  S,n"'(f)  .  0<  f-  <\ 

^  =  L  ±  ^  (?)> 

These  results  may  be  displayed  on  the  complex  plane  as  shown  in  Fig.  26. 


Figure  26.  Behavior  of  the  Characteristic  Exponent^x<«: 
as  a  Function  of  g* 


Thus,  stable  solutions  occur  only  for  ^  between  zero  and  unity.  The 
stability  of  the  moticn  may,  therefore,  be  partially  analyzed  without 
obtaining  an  explicit  solution,  merely  by  determining  the  conditions  under 
which  o  <  <5^  <  \  .  Unfortunately,  this  determination  requires  extensive 
numerical  calculations  since  5'^  depends  upon  ACO  anti  hence  upon  all  of 
the  Fourier  coefficients  . 

Before  discussing  these  numerical  calculations,  however,  it  is  instruc¬ 
tive  to  obtain  explicit  solutions  to  Eqs.  (83),  for  it  is  noted  that  once 
is  known,  the  solutions  for  ^  may  be  found  from  Eq.  (8li).  Confining 
attention  to  the  stable  case  vHnereln is  purely  imaginary  (call  it 
),  two  linearly  independent  solutions  of  Eq.  (83)  are  the  coiqplex 
functions 


y=!-oo 


(90) 
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Now  W  and  loj^  must  separately  satisfy  the  recurrence  formula  of 
Eq.  (86)  so  that, 

-  C (2w  ^ 


(91) 


fIS 


The  second  of  Eqs.  (91)  may  also  be  written  as 


00 

-  (2U  ° 


(92) 


«*-«o 


whereupon,  comparing  Eq.  (92)  and  the  first  of  Eqs.  (91),  one  concludes 
that  and  satisfy  the  same  recurrence  formula  and  hence. 


Th*  two  solutlont  given  by  Eqs.  (VC)  thcr  becoi 


=  X.  e 

*»=-•  (93) 

u^c<j3  =  ^  e 


The  unknown  { “*®}'  ^  found  from  the  firet  of  Eqa.  (91).  Since  this 
recurrence  foneile  ie  reel,  it  follows  that  the  [W,/}  ere  elao  reel. 
Conatent  ailtiplea  of  tho  eua  end  difference  of  the  two  aolutiona  ^.ven  by 
Eqe.  (93)  yield  respectively  the  two  linearly  independent  reel  eolutione 


@,C4>)2X(;u.+U*)=  ^ 

Jts-m 

end 

Pliullj^  th«  ooaplti#  hoaog«tnecu5  vlution  for  3  is  gisrn 


(94) 


§(<^)  =  ae,c<f>  T(^e*r4); 


where  CL  and  are  constants  (which  nuet  be  real  if  9  is  tr  be  a  reel 
function  of  the  reel  variable  ^  )  detemlned  by  the  initial  conditions 
of  the  Botion. 

Eqs.  (94)  show  that  the  solution  for  §  ia  wholly  periodic  of  period 
ir  or  2it'  only  when  is  an  integer,  which  Mane  that  these^  periodic 
solutions  occur  only  at  the  reapeetive  comers  ^  >  o  or  ^  •  C  of 
Pif.  26.  Thus,  the  unstable  rofimns  of  the  solution  ere  bounded  by  purely 
periodic  solutions.  In  the  ease  of  Hathieu's  equation  (where  only  a* 
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and  3i  differ  from  zero) ,  these  periodic  solutions  have  been  named  Mathieu 
functions  and  are  designated  as  even  or  odd  according  to  whether  they 
reduce  to  the  cosine  or  sine  function  as  a,— ►©  (Refs.  32-3U).  In  the 
general  Hill  problem,  the  corresponding  functions  do  not  appear  to  have 
been  tabulated. 

In  principle,  the  foregoing  analysis  yields  the  complete  homogeneous 
solution  for  &  ,  provided  that  the  exponent  may  be  determined  from 
Eq.  (87).  The  evaluation  of  which  is  required  for  the  calculation 

of  ,  however,  poses  certain  practical  problems  of  convergence  for 

moderate  to  large  values  of  Bo  ;  that  is,  for  large  values  of  .  The 

convergence  of  AC®!  may  be  stuped  the  analysis  which  follows. 

Consider  Eq.  (88)  for  Mathieu's  differential  equation  wherein  only 
So  and  Si  differ  from  zero.  Then  the  typical  term  in  the  Ptb  row 
of  /aCo)  adjacent  to  the  main  diagonal  (of  which  every  term  is  unity)  reads 


-ai 

-  a® 


all  other  terms  being  zero. 

Now  for  nximerical  purposes  the  infinite  Hill  determinant  is  replaced 
by  a  large  but  finite  determinant,  say  by  •  Let 

80 »  (.2. •ml'*'.  Near  the  center,  the  determinant  then  appears  as 

1  £0  o . 

e.  l£o . 

o  <£  I  e  . 


£  f  £  O 

o  £  I  e 

O  O  £  I 


since  is  negligible  compared  to  Bo  for  . 

A  determinant  with  such  symmetry  may  be  expanded  by  considering  the 
following  li  X  li  example: 
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J  £  o  o 
€.  c  £  o 


o  o  £,  \ 


Multiplying  the  first  row  by  S  and  subtracting  from  the  second. 


1  £  0  0 

1  £  0 

-£^  £  0 

0  i-e  £  0 

II 

0  €  1  € 

=• 

£  «  e 

+ 

0  1  6 

0  0  fi  1 

0  e  1 

0  £  1 

=  Dj  -  . 

This  relation  aay  be  readily  generalised  to  any  else  dstersdnant,  viz: 

D.  =  . 

Thus, 

D,  =1 

D^=  I  -z.6^ 

D4  = «  -  4-  6^ 

Dj  =  ‘  -  E*  ,  e^c.* 

which  may  be  further  generalized  to  read  (to  terms  of  order  ) 

D,  =  1  +  Q->'>e* 
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This  expression  shows  clearly  that  as  long  as  3o>'^  ,  examining 

successively  larger  sizes  of  determinants  id.ll  now  show  convergence  regard¬ 
less  of  the  size  of  <£  >  o  .  Convergence  does  not  begin  to  be  seen  until 
the  determinant  is  sufficiently  large  so  that  transcends  3©  , 

after  which  th*’  convergence  should  be  fairly  rapid.  The  result  also 
generalizes  to  the  Hill  problem  where  all  ■£3n^  may  be  present. 

Thus,  the  evaluation  of  Hill's  determinant,  which  is  troublesome  even 
with  the  aid  of  a  high  speed  computer,  prevents  the  wholly  analytical 
solution  of  even  the  homogeneous  problem  from  being  achieved.  For  small 
values  of  3©,  the  procedure  outlined  above  may  be  carried  out  manually. 
The  solution  of  the  non-homogeneous  equation  may  then  be  obtained  by  the 
method  of  variation  of  parameters  (Ref.  37,  Page  122).  The  technique 
should  yield  meaningful  results  if  applied  to  the  librations  of  the 
dumbbell-shaped  satellite  moving  on  an  eccentric  orbit,  but  for  the  case 
where  aerodynamic  influences  are  of  interest,  the  method  is  too  cumbersome 
to  be  of  direct  value. 

5 . 2  Stability  of  the  Natural  (Unforced)  Motion 

The  foregoing  analytic  solution  is  of  considerable  value,  however, 
in  providing  a  qualitative  description  of  the  stability  of  the  unforced 
angular  motion  on  an  elliptical  orbit.  Since  stable  motion  can  only 
exist  for  o  <  §  |  ,  it  follows  that 

O  <  S  (4Tr  )  <  1 


is  the  requirement  for  bounded  motion.  When  3*  is  not  the  square  of  an 
even  integer,  a  sufficient  condition  for  stable  (bounded)  motion  is 
o  <  <  I  }  however,  there  is  no  reason  to  expect  ^C«>>  to  remain 

on  this  interval. 

The  circumstances  under  which  the  above  inequality  is  violated  may 
be  examined  by  considering  the  integer  values  of  -vi  a.  .  In  cases  where 
a®  is  near  an  odd  square,  )  will  be  close  to  unity  vhile 

AC®>  will  be  determined  essentially  by  all  of  the  and  will  be 

relatively  insensitive  to  small  changes  in  a©  .  Regarding  A.CO  as 
independent  of  3©  for  the  latter  near  an  odd  squsure,  three  subcases 
result:  (i)  For  ACo^  <  I  ,  the  motion  remains  bounded  as  a©  approaches 
(zn+i)'*’  .  (ii)  For  A(o^  =  I  ,  a  single  unstable  point  occurs  for  each 
.  Since  I  in  this  case,  =  c  and  Eqs.  (91*)  show 

that  periodic  solutions  of  period  Z'Tr  exist,  (iii)  For  a  Co)  >  I  ,  a 
region  of  instability  may  be  associated  with  each  a®  »  .  At 

the  boundaries  of  this  region,  =  I  and  hence  periodic  solutions  of 
period  zTf  occur  on  these  boundaries. 
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The  case  of  ao*C'*.'oV*'  has  been  excluded  in  the  construction  of 
and  hence  this  case  requires  special  treatment.  It  is  not  difficult 
to  show  that  although  ^C«>)  is  undefined  at  ,  |p*approaches 

a  definite  limit  as  Unfortunately,  this  limit  is  not 

readily  calculable.  In  the  special  case  where  all  of  the  {3*,}  save  So 
vanish,  •=  |  and  a.  =  are  the  eigenvalues  of  the  sine 

problem.  Since  =  o  in  this  case,  Eqs.  (9ii)  show  that  periodic  solu¬ 
tions  of  period  it  result.  When  the  other  £dri}  are  present,  a  similar 
behavior  is  seen,  with  the  possibility  that  another  region  of  instability 
may  be  associated  with  each  3®= 

The  foregoing  discussion  brings  out  the  essential  features  of  the 
problem  of  stability  determination.  The  specific  behavior  may  be  described 
as  follows:  Let  all  of  the  {3»\]  be  fixed  with  the  exception  of  one,  and 
let  this  one  be  permitted  to  vary.  As  it  does,  will,  in  general, 

exhibit  an  oscillatory  behavior  which  will  cause  g'*-  to  leave  the  range 
of  o  <  g'*-<  I  .  Corresponding  to  ^^<c>  and  5*  >-\  ,  the  motion  will 
be  completely  unstable.  Corresponding  to  =  o  and  3*"=  1  ,  periodic 
(and  therefore  bounded)  solutions  of  period  tf  and  z-rr  respectively 
exist,  but  the  two  solutions  given  by  Eqs.  (9i:)  become  linearly  dependent 
in  each  case  and  the  second  independent  solution  will  be  found  to  be 
unbounded  (Ref.  32).  Thus,  alternate  regions  of  stable  and  unstable  motion 
will  generally  be  encountered,  with  each  unstable  region  bounded  by  periodic 
functions  of  period  Tr  or  z-rr . 

Since  this  behavior  is  evident  even  when  only  3®  and  3  |  differ 
from  zero,  it  is  convenient  to  continue  this  discussion  for  the  Mathieu 
equation.  The  stable  and  unstable  regions  for  this  problem  have  been 
described  by  numerous  authors  (Refs.  32-35)  and  may  be  displayed  on  a  plot 
of  So  versus  S|  (c.f..  Ref.  33»  Fig.  166).  For  the  present  case, 
however,  both  and  a,  (and  indeed  all  the  {.an"}  in  the  complete 
problem)  depend  upon  the  four  physical  parameters  cP^  ,  M  ,  e  and 
X  ,  as  seen  from  Eq.  (73).  For  illustrative  purposes,  assume  M  =  o  and 
neglect  g  in  comparison  to  CPp*  1,CX)  •  Then  Eq.  (83)  specializes  as 

e"  +  +  Z  ©  =  ■:- 

Thus ,  the  stable  and  \instable  regions  on  the  So  -  3 1  plane  may  be  mapped 
into  corresponding  regions  on  the  C?p*  -  X  plane.  Fig.  2?  presents  an 
approximate  construction  of  these  regions  based  upon  Ref.  33,  and  is 
intended  for  qualitative  examination  only. 

The  nature  of  the  stability  boundaries,  as  determined  from  Mathieu ’s 
problem,  appear  as  follows  for  a  given  configuration:  A  maximum  perigee 
altitude  exists  for  any  given  eccentricity  such  that  the  motion  may  be 


stable  everywhere^.  Within  this  stable  region,  however,  lie  isolated 
regions  of  instability.  In  some  cases  these  regions  reduce  to  merely 
lines,  particularly  for  large  values  of  (low  altitudes)  and  small 

values  of  X  (low  eccentricities),  but  in  other  cases  the  regions  may 
be  large.  A  typical  orbit  will  begin  at  some  value  of  hp  and  and 
€  will  start  to  decrease  while  hp  remains  fairly  constant.  During 
this  decay,  X  will  drop  rapidly  and  CPp^will  increase  (due  to  the 
increase  in  )  >  both  effects  tending  toward  greater  stability  of  the 
angular  motion.  However,  it  appears  quite  possible  that  these  changes 
will  carry  the  motion  into  a  small  but  finite  region  of  instability,  and 
that  the  vehicle  may  dwell  there  long  enough  to  begin  tumbling.  It  is 
also  possible  that  Xp  and  £  would  initially  place  the  vehicle  in  an 
unstable  region,  while  the  subsequent  orbit  decay  would  move  it  out  of 
this  region. 

These  results  correspond  to  intuition.  When  the  aerodynamic  stability 
margin  is  very  great,  the  motion  should  be  stable  as  it  is  at  low  altitudes. 
When  the  margin  is  small,  the  possibility  of  unstable  motion  exists, 
particularly  since  the  natural  frequencj*  of  the  aerodynamically  induced 
oscillation  (beii:g  of  order  >1  )  may  become  comparable  to  the  orbital 

frequency.  Furthermore,  for  a  given  value  of  the  stability  is 

greater  for  small  values  of  X  than  for  X  large,  as  might  be  expected. 

The  catastrophic  nature  of  the  amplitude  build-up  possible  in  these 
small  unstable  regions  is  illustrated  in  the  section  dealing  with  results 
of  the  numeric  integration  of  the  equation  of  motion,  and  hence  no  further 
elaboration  of  this  general  subject  appears  warranted  here.  It  is  suffi¬ 
cient  to  note  that  the  problem  of  passive  satellite  stability  appears  to 
depend  primarily  on  the  rate  at  which  the  amplitude  increases  when  an 
unstable  region  is  encountered. 

5.3  Stability  of  the  Forced  Motion 

The  preceding  section  has  shown  the  solution  of  the  homogeneous 
equation  to  contain  niunerous  unstable  regions.  The  addition  of  periodic 
excitation  (or  forcing)  will  not,  in  general,  mitigate  this  condition 
and  in  fact,  will  generally  compound  the  difficulties  by  permitting 
resonances  to  occur. 

An  analysis  of  the  forced  motion,  based  upon  the  properties  of  second 
order  differential  equations  with  periodic  coefficients  has  been  performed 
by  Juelich  in  Volume  III  of  this  series  of  reports  (Ref.  2)  and  need  not 
be  repeated  here.  The  essential  result  for  this  discussion,  is  that  the 
only  effect  of  the  periodic  excitation  is  to  induce  a  resonance  idien  the 


^For  M  >o  >  and  C?p'’*>o  f  maximum  altitude  lies  at  infinity 
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period  of  the  homogeneous  solution  and  that  of  the  forcing  function 
coin^idel.  In  terms  of  the  transformed  variable  <P  ,  this  occurs  when¬ 
ever  the  homogeneous  solu  bion  is  periodic  of  period  "O' ,  or  whenever 
=  o  .  Since  this  point  is  excluded  from  the  stable  range  of  the 
homogeneous  solutions,  no  additional  instabilities  in  the  motion  occur 
due  to  forcing. 

The  addition  of  periodic  forcing,  however,  does  permit  periodic 
solutions  of  period  If  (the  period  of  the  fbrcing  function  in  4  )  to 
be  found  for  all  non-resonant  conditions.  From  Ref.  2,  this  may  be  shown 
as  follows; 

Let  0o  be  that  solution  of  the  complete  equation,  Eq.  (73),  that 

has 


and  let  Of  and  be  those  solutions  of  the  homogeneous  equation, 
Eq.  (83)  that  have,  respectively. 


eiCo-i-  \  ,  ©/coi  «o 

so  ,  ©jf  Coi  rr  I 


(95) 


Comparing  these  initial  conditions  with  Eqs.  (9i;),  it  is  seen  that 


S’C'P)  = 

In  general,  any  arbitrary  solution  of  Eq.  (73)  may  be  written 

^C4>)  =  a  ©,C4>)  +  CB  -h 


(96) 


^hat  is,  when  the  system  is  "forced"  at  its  natural  frequency,  (the 
frequency  of  the  homogeneous  solution),  "resonance"  is  said  to  exist. 
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If  it  is  now  required  that  &  be  periodic  or  period  tr  (in  ^  ),  then 
it  may  be  shown  (Ref.  2)  that  suitable  constants  CL  and  CB  satisfy  the 
matrix  equation 


where  the  determinant  of  "W  is  the  Wronskian  of  the  solutions 
evaluated  at  <P  =  tr'  ,  and  I  is  the  identity  matrix.  A  unique  periodic 
solution  of  period  if  is  then  seen  to  exist  if  the  matrix  (l-W)  is 
non-singular.  With  the  help  of  Eqs.  (9U)  and  (96).  it  may  also  be  seen 
that 


and  so  (^I  -W )  is  singular  only  for  ®  >  which  has  been  identified 

as  the  resonant  case. 

The  existence  of  a  unique  periodic  particular  solution  for  all 
insures  the  existence  of  a  bounded  solution  of  the  complete 
problem  for  all  non-resonant  conditions.  In  the  ujistable  regions,  however, 
slight  deviations  from  this  periodic  solution  will  amplify  and  the  solu¬ 
tion  becomes  unbounded.  This  behavior  is  illustrated  in  Section  5.5.  It 
is  to  be  noted  that  as  resonance  is  approached,  the  periodic  solution  may 
grow  in  amplitude  so  that  in  the  neighborhood  of  a  resonance,  the  motion, 
although  bounded,  may  have  a  prohibitively  large  amplitude. 


Figure  27.  Stability  Boundaries  on  the  Op  -X 
Plane  Based  upon  Mathieu’s  Equation 
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S.k  Approximate  Solution  of  Complete  Pitching  Equation 

The  preceding  section  has  shown  that  in  general,  the  complete  solu¬ 
tion  of  the  equation  for  the  pitching  motion  of  an  aerodynamically  respon¬ 
sive  satellite  is  quite  complicated  and  extensive  numerical  computations 
are  required.  The  scope  of  these  calculations  is  such  that  a  solution 
may  be  obtained  with  little  more  effort  by  a  direct  numerical  integration 
of  the  equations  of  motion^  this  latter  approach  requiring  only  few 
simplifying  assumptions. 

It  is  of  considerable  value,  however,  to  have  available  a  completely 
analytical,  albeit  app:  >ximate,  solution  to  the  equations  of  motion.  Such 
a  solution  has  been  obtained  using  the  method  of  variation  of  parameters 
to  solve  a  reduced  version  of  Eq.  (73).  The  reduction  is  accomplished  by 
assuming  the  orbital  eccentricity  to  be  sufficiently  small  so  that__ 

and  all  ,  n  =  2.,  -  ■  multiplying  0  may 

be  neglected.  Thus,  the  homogeneous  equation  is  reduced  to  Mathieu's 
equation.  If,  in  addition,  all  terms  involving  except  the  first  are 

neglected  as  excitation  in  Eq.  (73),  the  equation  takes  the  form 

CP 

O  +  ^  c«-o-v^©  S  Dj,  S\n  r\  V"  (97) 

where,  from  Eqs.  (71)  and  (72), 

3  Ao  =  3  M  +  (P*  I„(\) 

^  =  (,I-3M>€ 

Z  Qp  Io(X) 

C,=  zf  [c?.*I,(X)-V] 

Eq.  (97)  is  in  the  standard  form  of  a  non-homogene ous  Mathieu  equation 
with  harmonic  excitation.  In  order  to  place  Eq.  (97)  in  the  customary 
form  far\a  variation  of  parameters  solution  (e.g..  Ref.  17),  it  is  neces¬ 
sary  to  duress  the  order  of  the  equation  by  writing 

=  ©  ,  02.=  (98) 


Then 


+  [a^  +  3-^'^]©,= 
0/ -  e^=o 


V4»l 


(99) 


are  the  pair  of  first  order  equations  equivalent  to  Eq.  (97). 

The  zeroth  order,  or  reference  solution  which  will  be  varied  is 
taken  as  the  circular  orbit  solution  of  Eq.  (77)? 

0=  a  Co-a  A  V  +<BSinA'xr  +  (lOO) 

where  the  "constants"  OL  and  <E»  will  be  varied  such  that  Eqs.  (99)  are 
satisfied  by  an  equation  of  tlw?  form,  of  Eq,  (100).  Thus, 

0^-0.  Cs.aA'vr  -+  <S>  S»n  A-ir  +  Pe> /A^ 


and 


Q  =  —dl  A  S»rv  A'U'  -h  (SA  AV“ 


0,  ^  =  d  ^  c«-o  Air  +C32>^SinAv*  -v  0^. 


(101) 


g)  '  =  _  <a|_'AS»n  A  V  -t  (B>^A  —  A^©>  ■+  P&  , 


Substituting  these  results  in  Eqs.  (99)  yields,  for  A  o  j 
CL  ^  C»-«  A  xr  +  ^^SvcAxT"  —o 

-CL*S\nAxr  +  <B'  Co-oAxr  =  -  i  Co.ovr^^^-^A  v-t  <S)Sia  AV^ 

CO 

4-  ^  Sin m-vr 

—  A  * 

n  =i 


These  equations  may  be  regarded  as  a  pair  of  non-htanogeneous,  linear, 
algebraic  equations  for  the  derivatives  d*  and  .  A  non-trivial 

solution  for  these  derivatives  is  assured  by  Cramer’s  rule  provided  that 
the  determinant  of  the  coefficients  does  not  vanish.  In  this  case  it  is 
clear  that 


93 


Co-oAv  SmAV- 


—  Si«Av  Co-o  Av 


I 


and  hence  a  solution  nay  be  found*  The  solutions  which  result  after 
suitable  algebraic  and  trigonometric  substitutions  are 

d'  -  SmC'2-A-Ovl  —  — ^  f  t«-oCe.A +0  "ir- 

•4  A  U  J  4  A  L 

OP 

+•  C»o(2.A-0‘'r  -  4  .^I^Cp-o  (.a  \r  -  Cp-d  C  A 


(10?) 


(^'  s-.-^L^Lfco^CZ-A+'O-U-  C^C*A-0V-  +  -Z-Cp^-U-l  --^SJ^S.eC^-A+'iir 
4  A  L,  .j  4  A  L 


4-  Si%%  C^A-0'y"J  +  z  CA  +  A’iir  -  Si^c  A-<^'> 


Eqs,  (102)  are  exact  at  this  point,  but  in  order  to  proceed,  it  is  neces¬ 
sary  to  perform,  the  integrations  to  obtain  ^Cv)  and  ©C'w)  .  Since  the 
quantities  a,  and  (B  appear  on  the  right-hand  sides  of  these  expressions, 
the  integration  must  involve  iteration.  If  G.  and  ©  do  not  vary 
rapidly,  however,  an  acceptable  first-ordor  approximation  may  be  obtained 
by  integrating  Eqs.  (102)  regarding  CL  and  ©  a»  constant  and  equal  to 
their  mean  values  oh  the  right  of  these  equations.  Thus,  for  AqA  o,  '/t  : 


a  -  a.,  =  - 


ZA^AA"*- 


C«w-oV'C<«ZAV  ~  \j  +  S<'%'U“  SintAV 


©  3-  r 

+  - 1  Z.A  Op-oir-  &nr*  t  AV“  4-  ^in"Vj"  (4A'^  —  I 

z.a(4a^-i)L  (103) 

y\  A  V  J 


4- 


,  n'vr  C’O-oAV" 


—  Y>  C.P-0  v»  -U"  Si» 


-■?! 

'I 

4 


9h 


CR-  ffi>„  =  —  — ^  ^ r 2  A  <>ov- +  SvA-vr  (4  A*”  —\  -  C«-o2,Axrbl 

®  2A<|A'--0L  J 


ca^  r. 


ZAC4A'--i)L 


-  2  A  (c-^ov- c-<r«  t  AXr  -  0  Sin  ir  Sno2AXr 

'•  J 


(lOA) 


CD 

+  /  - - -  AS‘'«v\-u'  Sin  A\r  4-Y\c*-ovw-  C^fiy\r  -  I  1 

^  ACA'--n'-)L  -1 


where  Cio  and  C&o  are  the  values  of  these  quantities  at  \r  =  o  and  where 
a.  and  <a  denote  mean  values.  Inserting  the  above  results  into  Eq,  (ICC) 
and  collecting  terms,  the  following  first-order  approximation  results; 


O  —  Ci - ^  ^  (co-0  v-l)  -¥■  Sin  ■u-l  AV" 

L  4A''-I^  4A'--l  J 


+  1©,.  -  Sin  nr  j A~u-  (105) 


+  — "  4-  V  - - -f  A  Sin  n■^J-  —  V\  Sin  A»xr~l 

A"^  Z_.  ACA^-n')L  J  • 

This  solution  is  completed  once  O.  and  Q6>  have  been  determined.  Applying 
the  mean  value  theorem, 


a 


I 

V 


■XT 


^  e-hc. 


N3  =  o,  1,  Zj 


to  Eqs.  (1C3)  and  (lOA)  with 


tr  •=  ■2.N  IK 


a[,- 


»  4VK1A  (A 


A*--')  -* 


do  - 


(3  ^  4A^->-> 

-iirwA.  (^a"*--* 


! —  r  j  — 


c,«^  ^^.  Aw-nr 


-  T'  — Ti- 

N  tr  ( A'- -«■*•)■*•!- 


C-*<  Z.AM'O' 


<gr|+— S: - ^ - -  SinAAN-rl 

L  4  A*- -I  4irNA  (^a'-- 0  -* 


-CB.  -■ 


r,  » 

MA  (4A»---l)‘‘ l- 


C^-o  Z.A 


Nir  J 


_  - r  HU  -  s.»>  ZAIUIV  1 


When  w»l  f  the?e  expressions  may  be  approximated  by 

H'-^]  =«» 

®['  *7^]  =®- 

Introducing  Eqs.  (1C6)  into  Eqs.  (105),  and  defining 


(1C6) 


s  =• 


a*  s. 


C2>*a 


4A*  -I 

0(e>>  -  Do /A*- 

I  -  % 


©Vo) 


ACl  +0 


(IC-?) 


% 


'here  results, 

OCv'i  ~  [^a*Cl-  StA-v-J  c^A- 


+  j^(£>^(l-  %Ctrair^  — '2Aa.*%  s»« 


Av 


(ICP) 


C0 
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where  second  order  terns  have  been  neglected  in  composing  Qo  .  The 
homogeneous  solutions  obtained  in  Eq.  (llO)  correspond  to  the  first  order 
solution  to  the  Mathieu  equation  obtained  by  Struble  and  Fletcher  (Ref.  35) 
using  a  modified  variation  of  parameters  approach.  In  addition,  for  the 
dumbbell  satellite,  the  solution  SoC'v^  becomes 

BoM  =  €[ —Sin'Ji’-v- (111) 

which  has  also  been  obtained  by  Frick  and  Garber  (Ref.  28). 

Eqs.  (110)  have  been  conpared  wilh  the  corresponding  solutions  deter¬ 
mined  by  direct  numeric  integration  of  the  "exact"  equation  of  the 
pitching  motion.  Two  cases  a^e  examined:  the  gravity  gradient  stabilized 
dumbbell  satellite  and  a  body  with  Psi  0.1  ,  M  =  I  and  =  o 
orbiting  at  200  mile  perigee  altitude.  The  results  are  presented  in 
Figs.  28  and  29  respectively,  for  small  eccentricity  orbits. 

Fig.  28  shows  that  the  approximate  analytical  results  accurately 
reproduce  the  numeric  cwiputations  for  the  dumbbell  satellite.  As  shown 
by  Frick  and  Garber  (Ref.  28),  the  solution  given  by  Eq.  (Ill) 

for  the  dumbbell  satellite  is  a  good  approximation  for  eccentricities  as 
high  as  one-tenth,  and  hence  the  results  presented  here  merely  confirm 
this  earlier  work.  When  aerodynamic  torques  are  present,  however,  the_ 
situation  is  complicated  by  the  appearance  of  an  infinite  series  in  • 
For  the  particular  conditions  used  for  Fig.  29,  satis^ctory  convergence 
was  obtained  by  using  only  the  first  three  terms  in  ^md  in  general, 

five  to  six  terms  should  be  adequate  for  even  the  largest  eccentricities 
of  interest.  Mote  from  Fig.  29(b)  that  although  and  are  dupli¬ 
cated  rather  well  by  the  approximate  solution,  shows  a  significant 

deviation  from  the  exact  results. 

In  general,  the  agreement  between  the  approximate  and  exact  results 
appears  sufficiently  good  so  that  Eqs.  (llO)  may  be  used  to  estimate  the 
effects  of  orbital  eccentricity  on  the  vehicle’s  librations,  at  least  for 
small  values  of  €'  .  It  must  be  noted,  however,  that  the  values 
A  =  o,  i/z.  ,  I  are  excluded  from  consideration  since  Eqs.  (110) 
become  singular  for  these  values.  Accordingly,  for  an  investigation  of 
the  natvire  of  the  motion  near  the  stability  boundaries,  numeric  integra¬ 
tion  appears  required. 
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5.5  Numerical  ReiBults 


5.5.1  General 


The  purpose  of  this  portion  of  the  study  is  to  determine  the  nature 
of  the  satellite's  pitching  motion  on  a  slightly  eccentric  orbit  idien  aero- 
cfynamic  and  gravitational  torques  act.  The  form  of  the  governing  differen¬ 
tial  eqi:iation  suggests  that  the  paramount  question  to  be  answered  is 
whether  or  not  the  regions  of  instability  suspected  to  be  imbedded  in  the 
region  of  satellite  motion  do,  in  fact,  exist  and  whether  or  not  these 
regions  actually  destabilize  the  vehicle.  In  Sections  5.2  and  5.3»  the 
analytic  solution  of  the  pitching  equation  is  discussed  and  a  stability 
dj.agram,  Fig.  27,  is  presented,  showing  that  for  \  small  and/or  (Pf* 
large,  stable  motion  is  to  be  expected,  \diile  for  X  large  and/or  (pp* 
small,  ims table  motion  will  generally  predominate. 

Inasmuch  as  this  analytical  treatment  relies  on  a  number  of  assump¬ 
tions  about  the  relative  magnitudes  of  the  various  factors  appearing  in 
the  complete  equation,  one  objective  of  the  numerical  analysis  is  to 
verify  the  existence  of  stability  boundaries  at  least  qualitatively  similar 
to  those  presented  in  Fig.  27.  Moreover,  the  numerical  analysis  is  also 
intended  to  substantiate  the  assumptions  and  order-of -magnitude  arguments 
made  in  arriving  at  the  approximate  equations  of  motion.  Finally,  th 
numerical  analysis  must  give  a  clear  picture  of  the  satellite’s  angular 
oscillations  as  a  function  of  time  and/or  position  on  the  orbit,  so  that 
frequencies,  amplitudes,  required  initial  conditions,  etc.,  may  be 
determined. 

It  may  be  asked  at  the  outset  if  it  is  possible  to  determine  the 
boundedness  of  the  motion  for  all  time  from  a  detailed  integration  of  the 
differential  equation  over  one  orbit.  The  results  of  an  analysis  of 
second-order,  non-homogeneous ,  linear  differential  equations  with  periodic 
coefficients  (Ref.  2)  indicate  that  such  a  determination  is  generally 
possible.  This  analysis  also  reveals  several  parameters  irtiich  may  be 
calculated  specifically  for  this  purpose. 

The  results  presented  in  this  section  are  based  upon  the  direct 
numeric  integration  of  the  "exact"  equation  of  pitching  motion,  determined 
in  Section  3.1.  Inasmuch  as  Ref.  2  contains  all  of  the  details  relative 
to  the  IBM  70li  computer  program  used  for  this  purpose,  only  the  essential 
elements  are  discussed  here. 

The  basic  differential  equation  which  is  solved  is  the  second  of 
Eqs.  (U6)  presented  in  Section  3.1.  This  equation  is 

i!f_  + x.csL+'fi’p  (112) 
dir^  dir  (.  J 


where 


2Q.=  2.  +€0-^^^:j]  -  2€Sm-Lr 

iQ  -  €  Smir 

°  I  +  (113) 


The  following  geophysical  constants  are  employed: 

1.U07,697,7  X  10^6  feetV  sec^ 

ae=  ^  “  20,902,131  feet(^) 

^3®t0g^»  7.292,115,1  X  10”^  radians/sec, 

and  the  atmospheric  density  variation  is  taken  from  the  1959  ARDC  (Minzner) 
atmosphere  (Ref.  26). 

The  computer  program  generates  the  three  solutions  to  the  exact  pitch¬ 
ing  equation  0,  ,  0,  and  ©•*.  (see  Section  5.3)  as  functions  of  xr- 

from  ir  =  o  to  In  addition,  the  constants  CL  and  (S  associated 

with  the  unique  periodic  solution  are  computed  for  those  cases  where 
Cx-W)  is  non-singular.  These  constants,  or  weights,  may  be  interpreted 
as  the  initial  conditions  which  are  required  in  order  that  the  motion  be 
precisely  periodic  and,  therefore,  bounded.  An  upper  bound  on  the  amplitude 
is  also  computed  for  each  orbit. 

The  basic  variables  vmdor  study  are  the  configuration  parameters  P  , 

,  M  and  the  orbit  parameters  (or  elements)  Vip  ,  e  and  6  .  For 
the  purposes  of  conqparison,  a  wide  range  of  these  parameters  are  investi¬ 
gated  in  various  combinations  using  the  digital  computer  program  described 
above.  A  small  analog  computer  mechanization  is  also  employed  for  check 
purposes.  The  results  of  these  studies  are  described  below. 


lou 


5.5.2  Diacussion  of  Reaults 

5. 5. 2.1  Range  of  Valuables 

The  six  significant  parameters  numerated  above  vere  studied  over 
the  following  ranges: 

o  <  p  <  \oo 

o  5.  £  1 000  ^+Vsla^ 

-I  £  M  <  1 

The  separate  effects  of  these  variables  are  discussed  in  the  following 
sections . 

5. 5. 2. 2  Effects  of  Aerodynaniic  Damping 

In  Section  U.3,  dealing  with  tlie  circular  orbit  results,  it  is 
shown  that  the  aerodynamic  damping  has  an  insignificant  effect  on  the 
vehicle’s  angular  oscillations.  In  the  current  numeric  investigation, 
values  of  up  to  100  fail  to  show  any  effect  idiatever  on  the  angular 
motion  on  the  circular  or  eccentric  orbits^.  Analog  computer  results 
likewise  reveal  no  appreciable  effect  of  aerodynamic  damping  at  altitudes 
as  low  as  100  miles  and  values  of  FJ  up  to  1,000. 

These  calculations  indicate  that  the  analytical  estimates  of  the 
time-to-danp  to  half-anplitude  (see  Appendix  C)  and  the  order  of  magnitude 
arguments  of  Section  3.3  are  substantially  correct.  Accordingly,  it  is 
concluded  that  the  aerodynamic  danq^ing  has  essentially  no  direct  effect  on 
the  vehicle's  angular  motion,  at  least  during  the  time  interval  for  ’irtiich 
the  orbit  elements  may  be  regarded  as  constants.  The  presence  of  even  a 
small  positive  damping  terra,  however,  indicates  that  the  motion  cannot 
grow  without  bound,  although  the  anplitude  may  become  excessively  large. 
For  the  purpose  of  examining  the  pitching  motion  further,  the  aerodynamic 
danping  will  be  taken  as  zero.  This  allows  the  vehicle  to  be  character¬ 
ized  by  the  two  configuration  parameters  P  and  M  , 

5. 5. 2. 3  Effects  of  Orbit  Inclination 


o  £  e  ^  O.OZS- 

XtS  ^  lip  ■s+ar+u+e  miles 

o  S  f  90“ 


The  effects  of  orbit  inclination  have  been  discussed  briefly  in 
Secti(xi  U.3.  Inasmuch  as  the  atmosphere  has  been  assumed  spherically 
symmetric  for  this  work,  the  entire  effect  of  i  is  taken  into  account  in 


^In  this  case  results  for 
decimal  places  printed. 


I 


»  O 


and  =  loo  agreed  through  the  four 
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the  term  Vr  /v  .  For  a  given  vsQ.tte  of  ^  ,  this  tern  may  be  regarded 
as  only  slightly  altering  the  values  of  P  and  Vj  and  accordingly,  no 
numeric  results  are  presented  for  this  portion  of  the  study. 

While  the  neglected  oblateness  of  the  avmosphere  produces  no 
major  direct  effect  on  the  motion,  it  does  Introduce  a  second  harmonic  in 
the  density  variationj  this  harmonic  having  period  TT  rather  than  period 
Ztr  .  Thus,  even  in  a  circular  orbit,  the  governing  equation  for  the 
pitching  motion  becomes  of  the  Hill  type,  and  in  view  of  the  analysis  of 
Sections  5.2  and  5.3,  unstable  conditions  may  exist  for  this  case.  This 
problem  appears  to  require  further  stucfy. 

5.5.2.U  Effects  of  Initial  Condlticais 

The  constants  CL  and  C£>  ,  determined  from  Section  5 .3,  may  be 
associated  directly  with  the  Initial  conditions  of  the  motion  and  Indicate 
the  importance  of  the  initial  conditions  in  composing  the  periodic  parti¬ 
cular  solution  of  the  equation  of  motion.  The  results  of  the  numeric 
calculations  show  the  solution  corresponding  to  0^  (that  solution  of  the 
homogeneous  equation  which  has  &,€«>:)=  i  ,  €#/Co>  =o)  to  be  given  an 
extremely  small  weight;  essentially  zero.  As  shown  in  Ref.  2,  ©p  must 
be  an  odd  function  whenever  Sc>  and  so  ©p  will  be  principally  odd 
for  small,  non-zero,  values  of  .  It  is  clear  that  if  the  periodic 

particular  solution  contains  ©»*with  a  non-vanishing  weight,  then  it 
cannot  be  odd.  Thus  0^  should  be  formally  absent  from  for 

5. 5. 2. 5  Effects  of  Remaining  Parameters 

The  parameters  P  and  M  influence  the  motion  in  essentially 
the  same  way  and  their  effects  cannot  be  separated.  Similarly,  the 
effects  of  Vip  and  €  are  most  conveniently  treated  together. 

Figs.  30  and  31  illustrate  the  behavior  of  the  pitching  motion  on 
the  elliptical  orbit  (G  =■  o.oisr>  at  UOO  miles  perigee  altitude  idien  P  or 
M  or  both  are  permitted  to  vary.  Far  Fig,  30,  only  a  very  small  amount 
of  aerodynamic  stability  is  present;  r=:C>.  lO  ,  Fig.  30(a)  shows  the 
solutions  for  M  =  o  and  Fig.  30(b)  shows  the  results  idien  Nd  —  I  .  In 
this  case,  the  significant  effect  of  M  on  the  amplitude  and  frequency  of 
the  oscillations  is  easily  seen.  For  Figs.  31(a)  and  (b),  on  the  other 
hand,  a  large  aerodynamic  stability  margin  is  present,  with  n=  loo  . 

Again  Fig.  31(a)  shows  the  results  for  M  =  o  irfiile  Fig.  31(b)  presents 
the  solutions  for  M  =  I  .  In  this  case  the  effect  of  M  is  substantially 
reduced.  This  general  behavior  characterizes  the  interplay  between  P 
and  M  throughout  the  altitude  and  eccentricity  range,  although  at  lower 
altitudes,  the  potential  effect  of  M  is  greatly  reduced  while  at  higher 
altitudes,  the  effects  of  even  large  values  of  P  may  be  readily  trans¬ 
cended  by  M  .  These  results  are  in  general  agreement  with  the  circular 
orbit  results  of  Section  It. 3. 
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Figs.  32(a)  through  32(g)  illustrate  the  typical  effect  of 
increasing  eccentricity  irfiile  holding  the  perigee  altitude  and  all  of  the 
vehicle  parameters  constant.  As  expected,  the  increased  eccentricity  both 
increases  the  amplitude  and  lengthens  the  period.  For  the  particular  set 
of  parameters  chosen  for  this  exanple,  an  unstable  condition  is  encountered 
between  €  =  o-ooz  and  e s 0.003  .  Further  discussion  of  these  unstable 
regions  is  presented  in  the  following  section. 

3. 5. 2.6  Stability  of  the  Motion 

The  numeric  integration  has  com'irmed  the  existence  of  stability 
boundaries  similar  to  those  presented  in  Pig.  27.  The  qualitative  behavior 
of  the  vehicle  as  it  passes  through  these  regions  of  stability  and  insta¬ 
bility  during  the  course  of  its  orbit  decay  is  described  in  Section  5-3 
and  need  not  be  repeated  here.  However,  it  is  of  interest  to  note  that 
the  length  of  time  the  vehicle  may  dwell  in  an  unstable  region  without 
catastrophic  amplitude  build-up  depends  upon  the  matrix  "W  of  Section  5.3 
and  upon  the  amount  of  deviation  from  the  periodic  particular  solution 
associated  with  the  orbital  conditions  (Ref.  2). 

In  order  to  illustrate  quantitatively  the  behavior  of  the  motion 
when  an  unstable  region  is  entered.  Figs.  33(a)  through  33(c)  have  been 
constructed  showing  three  complete  revolutions  of  a  satellite  on  its  orbit 
for  each  of  three  orbital  eccentricities,  all  other  orbit  and  vehicle 
parameters  being  held  constant.  In  each  figure,  the  dotted  curve  repre¬ 
sents  the  purely  periodic  particular  solution,  while  the  “solid^urve 
represents  the  motion  corresponding  to  a  2°  initial  misalignment  of  the 
vehicle.  The  three  eccentricities  could  well  be  encountered  by  a  satellite 
as  the  orbit  decays  due  to  aeroc^ynaitdc  drag. 

In  Fig.  33(a),  e.  =0.0  zs" ,  the  motion  is  stable  and  the  solution 
corresponding  to  the  initial  2°  misalignment  merely  oscillates  about  the 
periodic  particular  solution.  Fig.  33(b),  however,  is  the  same  vehicle  at 
an  eccentricity  of  e  =c>-oz.o,  which  lies  in  one  of  the  imstable  regions. 
It  is  noted  that  in  slightly  more  than  two  revolutions,  the  amplitude  has 
built  up  to  the  point  where  the  vehicle  will  certainly  begin  to  tumble. 

In  Fig.  33(c),  however,  the  unstable  region  has  been  passed  and  the  motion 
is  again  stable.  Upon  coitparing  the  solutions  for  the  first  half  revolu¬ 
tion  in  all  three  cases,  it  is  seen  that  they  appear  very  similar,  which 
might  lead  to  the  incorrect  conclusion  that  the  motion  for  €.=o.ozo  would 
be  substantially  the  same  as  for  e  =  «>.ozs'  or  e.-o.otS' .  Moreover,  the 
periodic  particular  solutions  appear  quite  similar  for  the  three  cases. 

Since  the  passive  oscillations  of  any  vehicle  permit  it  to  have 
an  infinity  of  possible  initial  conditions  upon  entry  into  an  unstable 
region,  it  is  clear  that  the  behavior  of  the  motion  shown  in  Figs.  33  must 
be  regarded  as  typical.  Thus,  the  situation  resolves  to  the  point  where 
passive  stabilization  of  artificial  satellites  depends  upon  the  rate  at 
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which  the  oscillations  grow  %hien  unstable  conditions  are  encountered  during 
the  orbit  decay.  These  unstable  regions  will  be  encountered  even  by- 
gravity  gradient  stabilized  satellites  of  non-spherical  shape  since  solar 
and  aerodynamic  eifects  produce  a  changing  center  of  pressure  which,  in 
general,  -will  not  coincide  with  the  vehicle's  mass  center.  The  situation 
can  be  improved  by  designing  the  vehicle  for  the  maximum  aerodynamic 
stability  possible  (see  Fig.  27),  but  the  isolated  unstable  regions  still 
exist,  and  must  be  considered. 

In  general,  the  design  of  a  passively  stabilized  satellite  will 
require  precise  knowledge  of  the  desired  orbit,  and  the  anticipated  orbit 
decay.  The  stabili-ty  boundaries  for  the  complete  problem  must  -then  be 
determined  In  a  form  similar  to  Fig.  27.  The  ex-temal  shape  and 

internal  mass  distribution  C  M )  of  -the  vehicle  may  -then  be  decided  upon 
such  that  the  vehicle  spends  its  useful  lifetime  in  a  stable  region.  An 
alternative  is  the  use  of  active  attitude  control  sys-tems  which  may 
shorten  the  vehicle's  useable  lifetime  or  compromise  some  of  its  other 
func-tions. 
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Figure  32(c) .  Behavior  of  Pitching  Motion  me  Orbital  Eccentricity 
Increases;  (c)  e^^o.oos' 
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Figure  33(a).  Pitching  Motion  on  an  Eccentric  Orbit  at  125  Miles  Perigee 
Altitude}  (a)  6=  o.ozs  ,  Stable  Motion 
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Figure  33(b).  Pitching  Motion  on  an  Eccentric  Orbit  at  125  Miles  Perigee 
Altitude;  (b)  €. -0.02.0  ,  Unstable  Motion 


figure  33(c) .  Pitching  Motion  on  an  Eccentric  l)rbit  at  12?  Miles  Perigee 
Altitude}  (c)  €.  =  0.0  1.5'  ,  Stable  Motion 
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SUMMARY  AND  CONCLUSIONS 


The  general  equations  governing  the  angular  motion  of  a  rigid 
satellite  are  derived  considering  the  orbital  motion  to  be  independent  of 
the  angular  motion  and  specified.  Both  gravitational  and  aerodynamic 
torques  are  included,  the  latter  being  developed  in  a  linearized  coeffi¬ 
cient  form  consistent  with  conventional  aerodynamic  approximations.  The 
atmosphere  is  assumed  to  be  spherically  symmetric  but  rotating.  Aerody¬ 
namic  damping  is  also  included  but  other  torques,  such  as  those  due  to 
solar  radiation  pressure,  are  considered  briefly  and  shown  to  be  negligible 
at  those  altitudes  where  aerodynamic  effects  are  significant. 

By  assuming  the  orbital  eccentricity  and  the  angular  deviations  from 
the  orbital  coordinates  to  be  small,  a  system  of  linearized  but  coupled 
equations  for  the  pitching,  yawing  and  rolling  motions  of  a  non-spinning, 
aid  symmetric  satellite  is  obtained.  These  equations  show  the  vehicle  to 
be  characterized  by  the  dimensionless  aerodynamic  parameters  SjT 

and  (which  describe  the  static  aerodynamic  restoring  moment  and  the 

aerodynamic  damping)  and  the  inertia  ratio  parameter,  M  (which  describes 
the  dyii«u'>ic  "shape"  of  the  vehicle) .  The  orbit  is  specified  by  the 
exponential  atmosphere  parameter  X  >  the  orbital  eccentricity  6  ,  and 
the  atmospheric  rotation  parameter  d  .  The  equation  for  the  vehicle's 
pitching  motion  is  shown  to  remain  independent  of  those  for  the  yawing  and 
rolling  motions  idien  aerodynamic  effects  are  considered.  The  rolling  and 
yawing  equations,  however,  remain  coupled  because  of  the  rotating  reference 
frame  employed.  The  angular  velocities  induced  by  pertiirbations,  such  as 
those  due  to  Earth  oblateness,  atmospheric  drag,  etc.,  are  found  to  be 
negligible  compared  to  the  angular  velocity  of  the  satellite's  motion  in 
the  orbit.  In  addition,  it  is  shown  that  for  small  eccentricity  orbits, 
the  ratio  of  the  velocity  relative  to  the  atmosphere  to  the  inertial 
velocity  remains  constant. 

The  resulting  equations  of  angular  motion  for  the  combined  influence 
of  aerodynamic  and  gravitational  torques  on  an  axisymmetric  satellite  are 
found  to  be  linear.  Order  of  magnitude  arguments  are  presented  which 
permit  the  equation  for  the  rolling  motion  to  be  integrated  directly  and 
the  equations  for  the  pitching  and  yawing  motions  to  be  greatly  simplified. 
A  Fourier  series  e:q}ansion  of  the  periodic  coefficients  of  these  equations, 
assuming  an  exponential  density  variation  and  small  orbital  eccentricities, 
yields  approximate  equations  of  the  non-homogene ous  Hill  type  which,  for 
very  small  eccentricities,  reduce  to  equations  of  the  non-homogeneous 
Mathieu  -type. 

The  equations  of  motion  are  solved  for  the  case  of  a  circular  orbit, 
and  both  qualitative  and  quantitative  descriptions  of  the  pitching,  rolling 
and  yawing  motions  are  found.  It  is  shown  that  in  this  case,  the  motion  is 
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con^letely  characterized  by  the  three  corfigxiration  variables,  C5p  , 
ol,*  and  M  ,  and  by  the  atmospheric  rotation  parameter,  ci  .  It  is 
shown  that  the  rotation  of  the  atmosphere  with  the  Earth  induces  a  steady- 
state  yawing  motion  of  the  vehicle  i^ich  prevents  yaw  equilibrium  from 
being  achieved  for  other  than  equatorial  orbits.  For  equatorial  orbits, 
a  completely  stable  equilibrium  is  found  to  be  possible  for  CPp*"  >  —  I 
or  (Pp*  >  3M  »  w]iichever  is  larger,  indicating  that  slender  bodies 
become  unstable  at  lower  altitudes  than  do  blunt  bodies.  Although  the 
pitching  and  yavdng  motions  may  be  rendered  stable  by  proper  design  of 
axisymmetric  vehicles,  the  rolling  motion  will,  in  all  but  exceptional 
circumstances,  conta.in  a  small  secular  term  which  will  cause  ^  to  even¬ 
tually  build  up  to  large  values.  This  effective  divergence  msy  be 
eliminated  by  relaxing  the  restriction  to  axisymmetric  shapes  or  by 
providing  active  roll-yaw  stabilization. 

The  effects  of  the  various  configuration  and  orbit  variables  are 
studied  for  the  circular  orbit  case.  It  is  found  that  gravity  gradient 
effects  will  generally  not  dominate  the  motion  until  altitudes  in  excess 
of  30O-UOO  miles  are  achieved,  while  at  lower  altitudes,  even  small 
amounts  of  aerodynamic  stability  will  suffice  to  override  the  gravitational 
toixiues.  The  aerodynamic  damping  and  the  orbit  inclinatJon  are  shown  to 
have  only  secondary  effects  on  the  pitching  motion. 

The  general  pitching  motion  on  an  eccentric  orbit  is  investigated 
based  upon  the  properties  of  Hill's  differential  equation.  The  exact 
solution  is  discussed  and  an  analysis  of  the  stability  of  the  solution  is 
performed.  It  is  foxmd  that  numerous  unstable  regions  may  be  encountered, 
in  which  any  small  deviations  from  a  purely  periodic  motion  (of  period 
equal  to  the  orbital  period)  will  amplify.  Thus,  stable  motion  exists 
only  for  certain  specified  combinations  and  values  of  the  parameters  P  , 
M  ,  and  €  .  The  periodic  forcing  or  excitation,  which  is  a 

natural  consequence  of  combining  aerodynamic  and  gravitational  torques,  is 
shown  to  permit  a  resonance  phenomenal  to  occur  on  the  boundaries  of  the 
previously  found  unstable  regions. 

The  results  of  this  analysis  indicate  that  the  ability  to  passively 
stabilize  an  artificial  satellite  rests  on  whether  or  not  the  orbit  decay, 
produced  by  atmospheric  drag,  causes  the  vehicle  to  enter  one  of  the 
unstable  regions  and  to  dwell  there  sufficiently  long  for  the  angular 
oscillations  to  grow  excessively  large.  Some  results  of  extensive  numeric 
integrations  of  the  equation  of  pitching  motion  are  presented  to  illustrate 
the  catastrophic  nature  of  these  unstable  regions  which  may  appear  in  a 
situation  giving  the  outward  appearance  of  stability.  In  this  connection, 
the  sensitivity  of  the  motion  to  the  initial  conditions  is  studied.  An 
approximate  solution  for  the  pitching  motion  on  an  elliptical  orbit  is  also 
obtained  using  the  method  of  variation  of  parameters.  The  results  are 
conpared  with  a  numerical  integration  of  the  complete  governing  equation. 
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It  is  concluded  that  the  aerodynamic  influences  play  a  decisive  role 
in  the  angular  motion  of  an  Earth  satellite,  whether  or  not  the  vehicle 
is  intended  to  be  aerodynamically  stabilized.  Since  it  will  be  virtually 
impossible  to  insure  that  the  center  of  pressure  and  the  mass  center 
exactly  coincide  on  any  given  satellite,  even  these  small  residual  aero¬ 
dynamic  torques  constitute  what  is  potentially  a  major  disturbing 
influence  on  the  satellite’s  orientation.  Through  proper  vehicle  design, 
however,  these  aerodynamic  influences  may  be  used  to  achieve  partial  or 
complete  attitude  stabilization  and  control.  In  particular,  the  design  of 
a  satellite  for  maximum  aerodynamic  stability  appears  to  reduce  the 
likelihood  of  encoxintering  an  unstable  region,  at  least  for  low  eccentri¬ 
city  orbits. 

The  need  for  additional  and  immediate  effort  in  this  general  area  is 
indicated  by  these  results.  Two  problems  of  significance  are:  (l)  the 
further  study  of  the  properties  of  the  non-homogeneous  Hill  equation  found 
to  govern  the  satellite's  pitching  motion  (in  order  to  obtain  the  stability 
boundaries  with  a  high  degree  of  accuracy),  and  (2)  the  synthesis  of  active 
aerodynamic  and  inertial  attitude  control  systems  based  upon  the  results 
presented  herein. 
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APPENDIX  A 


GRAVITATIONAL  TORQUES  ON  A  RIGID  BODY 


Consider  a  collection  of  naas  points  yf\i  viewed  froo  an  inertial 
refere^nce  frane.  Denote  the  Inertial  position  vector  to  the  1th  particle 
by  .  Then  the  equation  of  notion  for  this  1th  particle  Is  simply 


F"-  =  YVl  •  Q . 


(A-1) 


Now  If  the  1th  position  vector  Is  cross  multiplied  Into  each  side  of  this 
equation,  and  sum  taken  over  the  n  particles  of  the  collection,  the 
result 


f  5;  X  ^  =E  3,. ?,  =  i  f 


At  4- 


(A-2) 


Is  obtained.  Next  define  the  mass  center  of  the  collection  In  the  usual 
way  (Ref.  31)  such  that 


Z  1;  "i  =° 


(A-3) 


where  2^-  Is  the  position  vector  of  the  ith  particle  measured  in  axes 
fixed  in  the  mass  center.  Denote  the  inertial  position  vector  of  the  mass 
center  by  r  .  Then  clearly 


Q,  =  r  +  li 


and  the  left-hand  side  of  Eq.  (A-2),  which  is  the  net  torque  on  the 
collection,  becomes 


IS/ 


Also,  since  the  torque  due  to  ^e  internal 

forces  cancel  in  pairs  and  the  first  sunmation  Is  simply  'r  x 
the  moment  of  the  total  external  force  about  the  Inertial  origin.  'Hie 
right-hand  side  of  the  equation  of  motion,  Eq.  (A-2),  expands  as 
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^  E  +1i) 


t»‘ 


where 


•+  r)«2]7- Virii  +  Mr^r’ 1 

i*\  *  J 

Kl 

m=2I 


(A-4) 


(A-5) 


and  all  time  derivatives  of  are  relative  to  inertial  space.  By  the 
well-known  Coriolis*  law  (Ref,  31 )»  these  inertial  time  derivatives  may  be 
rewritten  in  terms  of  derivatives  with  respect  to  the  rotating  frame  fixed 
in  the  collection  with  origin  in  the  mass  center,  as 


2;  SI*  k 


so  that 


(A-6) 


For  a  rigid  collection  -o  and  by  Eq.  (A-3)  the  last  term  in 
Eq.  (A-6)  vanishes.  For  the  same  reason,  the  second  term  on  the  right  of 
Eq.  (A-4)  also  vanishes.  Finally,  the  last  tern  on  the  right  of  Eq.  (A-4) 
is  simply 


since  the  total  mass  of  the  collection  is  assumed  to  be  constant.  The 
coiqplete  equation  of  angular  motion  is,  therefore. 
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"!•  r 


(A-7) 


The  first  and  last  terms  in  this  expr^sion  cancel  since  together  they 
simply  represent  the  position  vector  r  ,  cross  multiplied  into  Newton*s 
second  law  written  for  the  mass  center  of  the  collection.  The  only 
remaining  term  on  the  right  of  Eq.  (A-7)  is  defined  to  be  the  true  rate  of 
change  of  angular  momentum  of  the  collection  viewed  from  the  moving  axes: 


Expansion  of  this  term  yields  the  acceleration  side  of  the  well-known 
Euler  equations  of  angular  motion  for  a  rigid  collection  (Ref.  31)  in 
terms  of  quantities  expressed  in  axes  fixed  in  the  collect.,  on.  The  princi¬ 
pal  problem  which  remains  is  to  expand  the  external  torques  for  a  rigid 
collection  moving  in  an  essentially  inverse-square  attracting  force  field. 
The  external  force  due  to  gravitational  attraction  on  any  mass  particle  in 
the  collection  is  simply 

pCs^T)  __  vifj.  I 

L  ^ 


where  §  is  the  potential  function  for  the  force  field.  The  torques  due 
to  forces  internal  to  the  collection  have  already  been  observed  to  cancel 
in  pairs  and  hence  only  the  external  forces  are  of  concern.  For  an 
inverse-square  attracting  force  field. 


where  the  origin  of  the  inertial  axes  are  here  assumed  to  coincide  with 
the  force  center  and  is  the  product  of  the  universal  gravitational 
constant  and  the  mass  of  the  primary  body.  Thus,  the  gravitational  force 
on  the  1th  particle  is  directed  along  the  inertial  position  vector,  : 
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and  hence  the  vector  torque  due  to  gravity  about  the  maM  center  of  the 
collection,  which  is  denoted  by  ,  is 


or 


N 


4*» 


=  ^2 

/  Bl  I 


=  -^11 


li 


(A.8) 


Clearly  the  last  tera  vanishes  from  the  definition  of  the  vector  product 
and  the  remaining  term  on  the  right  of  Eq.  (A>8)  becomes 


n 


2^  yfiL 


(A-9) 


The  suaaatlon  does  not  vanish  by  the  center  of  mass  condition,  since 
may  be  interpreted  as  a  weighting  factor  which,  in  general,  does  not  permit 
cancellation  of  the  torques  due  to  synaetrlcally  placed  mass  elements. 

Now  by  definition. 


=  r^+  l?^  +  zr.t 
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and  so 


3/z 


If  the  maximum  dimension  of  the  collection  i*'  small  compared  to  r*  ^ 

is  small,  and  i/o?  may  be  expanded  by  the  binomial  theorem  as, 


—  =  — (l-3— ‘ 
r3l 


and  hence  Eq.  (A-9)  becomes 


Again  the  first  sunaiation  vanishes  by  the  center  of  mass  condition, 

Eq,  (A-3),  while  terms  of  order  (Z;/r)^  may  be  neglected.  For  "small” 
bodies,  then,  the  torque  due  to  the  inverse-square  gravitational  attrac¬ 
tion  is  simply 


3 


^  r 


Z 


(A-IO) 


This  result  has  also  been  recently  obtained  by  Nidey  (Ref.  20), 

Up  to  this  point,  nothing  has  been  said  about  orientation  angle  and 
hence  no  small  angle  restrictions  are  inherent  in  this  result.  The  ^ 
detailed  nature  of  Eq.  (A- 10)  may  be  further  examined  by  expanding  r 
and  Ci  in  terms  of  components  in  the  moving  axes.  Asstime  rectangular 
axes  fixed  in  the  mass  center  of  the  body.  Then  let 


^For  a  collection  of  maximum  extent  10^  feet  at  the  Earth*s  surface,  Z^/r  ~  lo. 


(A-11) 


r=  2 


and  I. 


80  that 


9 

=z 


and 


+  (A-12) 


Inserting  Eqs*  (A-ll)  and  (A-12)  into  Eq»  (A-10),  the  gravitational  torque 
components  about  the  body  axes  become 


^3^.  " 


-  3  TT-Z  C^>  ^  Wi 

t=' 

-t  **  * 

t-i 

►L, 


Expanding 


t»‘ 

-nn=CiX;^  -^n=Cil  -ti'Xi.Ii,}  W; 
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Defining  the  moments  and  products  of  inertia  in  the  usual  way  by 


and 


_  N 


then  the  gravitational  torque  component  about  the  tc,  -body  axis  is 


'■*  '■>D  +''-zr,(l^-l,)}  (A- 


13) 


In  a  similar  way  the  other  two  gravitational  torque  components  are  found 
to  be 


+r;Yi(Vljj 


(A-U) 


!f' 

It  is  generally  most  convenient  to  let  the  body-fixed  axes  be  so  chosen 
that  the  products  of  inertia  vanish  (i.e.,  principal  axes)  in  which  case 
the  torques  assume  the  simple  matrix  form 


(A-15) 


Note  that  there  is  no  torque  about  any  axis  of  symmetry.  The  terms  ri/r~ 
are  merely  the  direction  cosines  between  and  the  instantaneous  body 
axes,  and  may  be  small  or  large  as  required.  Eq.  (A-15)  readily  reduces 
to  the  results  of  Roberson  (Refs.  9,  10)  when  the  Euler  angles  used  in 
that  reference  are  employed. 
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APPENDIX  B 

EXPANSION  OF  COEFFICIEKTS  IN  EQUATIONS  OF  ANGULAR  MOriON 


It  Is  well-known  that  the  unperturbed  motion  of  a  particle  in  an 
InTerae-equare  central  force  field  remains  planar  and  occurs  along  one  of 
the  conic  sections  with  the  force  center  at  one  focus.  The  polar  fonn  of 
the  equation  may  be  written  as 


where  r  is  the  distance  from  the  force  center  to  the  particle,  ~f  is 
the  semi-parameter  or  semi-latus  rectum  of  the  path,  £  is  tlie  eccen¬ 
tricity  and  nj-  is  the  true  anomaly.  Fig.  B-1  should  make  the  geometry 
clear. 


Figure  B-1 

The  form  of  the  path  depends  upon  the  value  of  ^  according  to  the 
following  table: 

|€'|>(  ,  the  path  is  a  hyperbola 
S'!  >  *  parabola 
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}6l<i(  j  the  path  is  an  ellipse 
\€)  =  o,  the  path  Is  a  circle  . 

For  the  present  purpose,  only  closed  paths  (hereafter  called  orbits)  are 
of  interest  and  so  the  restriction  nay  be  inposed* 

From  the  sketch,  the  semi-major  axis  (or  mean  distance)  of  the  ellipse 
has  length 

a  =  i  (rj  +  rp) 


and  so  it  nay  be  seen  that 


T  =  aC>-  6^;  . 


The  eccentricity  £  and  the  mean  distance  a  are  usually  the  fundamental 
constants  of  the  orbit.  In  the  present  case,  however,  it  is  convenient  to 
rewrite  Eq,  (B-l)  for  closed  paths  as 


r  = 


rp  ( I 

I  -hS  v 


\e\c\  . 


(B.2) 


Eq.  (B-2)  completely  defines  the  orbit  size  and  shape. 

The  velocity  vector  of  the  particle  is,  by  definition,  always  tangent 
to  the  orbit.  Defining  the  angle  between  the  normal  to  the  radius  vector 
and  the  velocity  vector  as  -y*  (see  sketch)  then 


r  v  =  V  Oo-o  "r* 


f"  =  V  Sio  V* 


and 


T a«  IT*  — 


r 

r^o- 


(B-3) 


Differentiating  Eq.  (B-2)  with  respect  to  time. 


13b 


•  te  Y"p  C I  • 

r  =  - - -r  VSirtir 

(\  + 


or 


•  I  •  £ 

r  =  V - Sio  ^r 


rpCi+-6> 

and,  combining  Eqs,  (B-3)  and  (B-4)  leads  to 


G  Sinxr 
\  *+•  S  C^oV 


which  Implies 


^\r^yr*  - 


<4-  _  G  S»n  v~ 


"/\  ^  «1.S  C©-a  ~ij~ 


c<^r*= 


TM  '+2.6  -v"  -V-  g'' 


and 


V 


=  n  Z.& 

Tp  V  0 +€'>*■ 


<l4ro-U" 


(B-4) 


(B-5) 


Making  use  of  the  fact  that  motion  in  an  inverse  square  force  field 
conserves  angular  monentum  relative  to  the  Inertial  coordinate  origin, 

r'*'v'  *»  H  p  a  c^«5+av\-V  (B-6) 


or. 


r"*-  V  =5  r-V  6o-o  TP*  =  XT 

p 


(B-7) 


and  so 
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^  _  jr  _J _ 

Frcu  Eqa.  (B>2)  and  (B-5): 

V  _  0+«)  ■/i'  +  26Co«'vr  +€'*■ 

Tp-vr  (^l  ■+€C-o^'\r)^ 

Again,  using  Eqa,  (B-2)  and  (B-7) 


V 

rp^ir 


Differentiating  with  respect  to  tiae 


or, 


V  =  2. 


\  4'ec<>t»xr 


0  +-6) ' 


(j-  6  Si«  v)  XT  v-j. 


•  • 

XT 


6  Sin  Xr 

z - 

I  4-feCi»-o-\j 


(B-8) 


(B-9) 


(B-10) 


The  coefficients  in  the  equations  of  motion  developed  in  the  text  may 
now  be  expanded  in  terms  of  xr  ,  the  independent  variable  and  and 
€  ,  the  constants  of  the  orbit*  Using  the  definitions  of  <Pp  ,  <Slp  ^nd 
(R,  as  given  by  Eqs.  (45)  and  the  definition  of  P  ,  and  r  as 

given  by’ £qs.  (46),  it  is  trivial  to  find 


ZR  =  a<,,(^)(^) 


V 

t>  XT 


V 


r\><r 


XT 


xr 


XT 


XT’ 


(B-11) 


Combining  Eas«  (B-8)  and  (B-IO)  with  Eqs.  (B-ll)  yields  the  first  three 
of  Eqs.  (47;  in  the  text. 

The  small  angles  A  and  TT  are  given  respectively  by  Eqs,  (36)  and 
(37)  in  the  text.  These  may  be  rewritten  for  small  values  of  £.  with  the 
help  of  Eq,  (B-7)  as 


S « 'f'  L  C.-0-0 


U. 


and 


v-p  \  V 


Introducing  Eqs.  (B-9) 


coa 


"vr 


S\rt  L  Co-oix  ] 


l+G 


(.  I  -4-  C  ) 


and,  with  the  aid  of  Eqs.  (B>5)» 


\  +  GCoo-u-  (.  "U-p.  I_l  -h  6  t-o-ovj  j 


(B-12) 


(B-13) 


As  the  final  development  in  this  Appendix,  consider  the  term 
^/r5-u-*.  From  the  definition  of  f  ,  and  Eqs,  (B-l)  and  (B-6) 


T 


^4 


and  hence 
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APPENDIX  C 

EFFECTS  OF  AERODYNAMIC  DRAG  AND  DAMPING  ON  THE  MOTION  OF 
NEAR-EARTH  SATELLITES 


1.  Drag  Effects 

A  primary  problem  in  the  study  of  satellite  aerodynamics  is  the  effect 
of  perturbation  forces  on  the  size  and  shape  of  the  orbit,.  Since  atmos¬ 
pheric  density  is  a  dominant  influence  on  aerodynamic  stabilization,  these 
changes  in  satellite  altitude  (and  hence  density)  are  very  significant; 
consequently,  the  effect  of  drag  perturbations  on  orbital  eccentricity  and 
altitude  decay  must  be  given  initial  consideration.  The  approach  to  this 
problem  is  to  assume,  as  a  first  approximation,  that  only  the  atmospheric 
retardation  and  the  attraction  of  a  spherical,  homogeneous  Earth  act  on  the 
satellite.  Furthermore,  the  drag  of  the  satellite  is  assvuned  to  vary  only 
with  air  density  and  the  square  of  the  relative  velocity;  1.  e.,  the  drag 
coefficient  is  constant  throughout  the  orbit. 

The  perturbations  which  the  drag  force  produce  on  the  orbit  elements 
may  be  determined  from  the  classical  method  of  variation  of  parameters 
(see,  for  example,  Moulton,  Ref.  17).  Sterne  (Ref.  19)  has  applied  this 
theory  to  determine  the  secular  variations  of  orbit  size,  shape  and  spatial 
orientation  for  a  satellite  with  drag  moving  through  a  rotating  atmosphere. 
For  the  present  case  a  static  atmosphere  will  suffice  and  Sterne’s 
equations  reduce  to 


—  =  ctE 


AL  =  -ib 

it  it 


(C-1) 


where  a  is  the  mean  distance,  €  the  orbital  eccentricity,  E  the 
eccentric  anomaly  and  b  the  ballistic  coefficient.  The  relevant  geometry 
is  shown  in  Fig.  C-1.  In  order  to  perform  the  integration  of  Eqs.  (C-l), 
it  is  necessary  to  express  the  density  in  terms  of  3  ,  6  and  E  .  The 
exponential  representation  given  in  Eq,  (52)  of  the  text  provides  the 
required  relation.  Using  the  eccentric  anomaly  in  order  to  rewrite  the 
equation  of  motion,  there  results. 


p-  =  a  -  6  co-o  E  ^ 
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and  so,  Eq,  (54)  of  Section  3.2  becomes,  without  further  approximation, 

^  ^  e  (c-2) 

where  X=  S  and  is  the  density  at  r*=  a  (or  E  «=  'TT /2_  ), 

Siuice  the  present  analysis  considers  only  orbits  of  small  eccentricity 
(say  less  than  O.l),  Eqs.  (C-l)  may  be  simplified  with  little  loss  in 
accuracy.  If  terras  of  order  and  smaller  are  neglected,  then  Eqs.  (C-l) 
and  (C-2)  may  be  written  as 

da  =  Z6c<K,E)ciE 

(C-3) 

d€  =  -Zba  (l-€^)  +6C^e)o>oE  dE  , 

The  changes  in  a  and  e  may.  now  be  found  by  integrating  Eqs.  (C-3) 
over  E  ,  considering  the  orbit  elements  appearing  on  the  right-hand  sides 
as  constants.  The  integrands  may  be  expanded  in  Fourier  series  using  the 
technique  discussed  at  length  in  Section  3.4.  As  before,  frequent  use  must 
be  made  of  the  integral  form  of  the  modified  Bessel  functions.  Performing 
the  integration  from  E  =  E«  to  any  E  and  confining  attention  to  the 
secular  terms  (those  terms  which  are  directly  proportional  to  E  ),  there 
is  obtained. 


e-6,=  -2ta0-e'>i°,[ux)  *6l.W  -f  I,W](E-E.)  (c-t) 

-i"  ^Pet-ioAic.  'T"eT-»r\S^  , 

Since  a  satellite’s  altitude  Is  minimum  at  perigee,  it  is  usually 
more  convenient  to  use  the  change  in  perigee  radius  rather  than  the  change 
in  send. -major  axis.  This  conversion  may  be  obtained  from  the  geometric 
relation  ir^=  aC'~6)  which  yields 

Arp  ^  -  a^€  .  (C_5) 
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f 


4 

f 


By  substituting  Eqs.  (C-4)  lnt.o  Eq.  (C>5)  and  neglecting  terms  of  the  order 
e/X  and  smaller^ 

^  (i-ie:>[r/V)-I,(\)]CE-EO  .  (C-6) 


The  second  of  Eqs,  (C-4)  and  E  q.  (C-6)  may  be  written  in  terms  of  the 
perigee  radius  as, 


>-p  -  rp.  3  -  2.  b  -  1,IV)]  (,£-£„> 


(C-7) 


6-e.  s  -  2.b  ^pO+«)l“^Lr,U)  7  6 1,Wl(E-E.) 


where  (|-ze)  =  0-€)^  and  where  (£/X)T^(,^)  is  neglected. 

The  only  remaining  problenmin  using  Eqs,  (C-7)  is  the  determination  of 
the  modified  Bessel  functions  I.  auid  X,  .  Since  the  argument  X  is  a 
function  of  altitude  and  since  1,  and  I,  vary  essentially  exponentially 
with  'h  ,  the  magnitude  of  and  (r»  -  are  highly  sensitive 

to  the  selection  of  a  references  altitude#  Part  of  the  problem  may  be 
alleviated  by  noting  that  the  mverage  density  over  m  complete  revolutions 
is 


zNtr 

ENir  r*> 

WA 


or 

1.0=  F/p,  • 


Now  if  the  average  density  is  aatumed,  as  a  first  approximation,  to  be  the 
arithmetic  mean  of  the  perigee  wtd  apogee  densities,  then 


zr.w  = 


Pe 


(C-8) 


The  advantage  of  Eq.  (C-8)  is  that  Io(V^  niay  be  conqsuted  without  know¬ 
ledge  of  an  effective  \  and  since  I,  is  a  known  function  of  I  „  > 

Eqs.  (C-7)  may  be  solved  (for  a  given  ballistic  coefficient)  in  terms  of 
initial  orbit  geometry  and  densities.  It  should  be  noted  that  for  eccen¬ 
tric  orbits,  the  density  at  the  average  altitude,  ,  is  always  less 

than  the  arithmetic  average  density,  i.e.,  initially 

circular  orbits,  X«C'°)  =  1.0  and  I,£.<5>=  o  which  leads  to  considerable 
simplification  of  the  first  of  Eqs.  (C-7)  along  with  the  deletion  of  the 
second  of  Eqs.  (C-7). 

Tnis  brief  treatment  of  drag  effects  on  satellite  orbit  decay  shows 
that  variations  in  the  ballistic  coefficient  may  alter  the  results  consid¬ 
erably,  since  the  decay  of  both  altitude  and  eccentricity  are  directly 
proportional  to  this  qxiantity.  The  problem  in  estimating  b  is  not  neces¬ 
sarily  attributed  to  the  drag  coefficient  alone,  but  to  the  effective  drag 
area.  Sn-  ,  upon  which  Co  is  based.  Calculations  readily  show  that  the 
lifetime  of  a  high  fineness  ratio  satellite  is  appreciably  reduced  if  it 
tumbles  rather  than  remaining  in  a  streamlined,  stabilized  condition. 


Effects  of  Darning  in  Pitch 


In  estimating  the  effects  of  pitch  damping  on  the  satellite's  oscil¬ 
latory  motion,  it  is  convenient  to  examine  the  time  to  damp  an  initial 
oscillation  to  one-half  its  initial  amplitude.  From  Eqs.  (62)  this  decay 


is  simply  expressed  as 


-  \  QC-fciJt 
0  J-o;  =• 


=  '/z 


(C-9) 


where  t'e/j.  denotes  the  value  of  the  true  anomaly  for  which  the  oscilla¬ 
tion  amplitude  has  been  halved  and  where 


—  e-sinv- 

p  (irp  (ytt\ 


(C-IO) 


Substituting  the  first  of  Eqs.  (C-IO)  and  taking  logarithms  of  both  sides. 


2p  jj  -  esin-fc|cl-t  =  0.693  . 


(c-11) 


This  integration  may  be  performed  by  following  the  procedure  outlined  in 
the  preceding  section.  The  result  is 


Sp*"  -  0.69>3 


(C-12) 


where  C^/t  ~  is  proportional  to  the  number  of  revolutions  required 

to  damp  to  one-half  the  initial  amplitude.  Initial  values  are  assigned  to 
e  ,  ^  and  a*  as  approxi.mations  since  orbital  decay  affects  these 
parameters  between  the  limits  of  integration. 

In  order  for  the  damping  parameter  to  have  a  significant  influence  on 
the  decay  of  the  angular  motions,  the  time  to  damp  to  half  amplitude  must 
be  considerably  less  than  the  time  required  for  drag  to  produce  a  signifi¬ 
cant  change  in  the  orbit  elements  themselves.  A  relative  measure  of  these 
two  times  may  be  found  in  the  ratio  of  from  Eq.  (C-12)  and 

the  required  for  a  change  from  Eq,  (C-7).  That  is, 


"^/zr 


Nez-g.  ^  b 


(C-13) 


where  N^v-  is  the  number  of  revolutions  taken  in  descending  feet 

(a  positive  increment),  and  where,  to  terms  of  order  £  ,  the  density  at 
r  =.  a ,  Pte ,  is  interchangeable  with  the  density  at  the  semi-latus 
rectum,  Pji  With  the  definition  of  the  ballistic  coefficient  and  the 
second  of  Eqs,  (C-IO),  Eq.  (C-13)  becomes 


=  11.09 


l.W  -  l,tV)  ]  r. 


(^C' •*■€:)  Io(X)  “  I^C^) } 


or,  conservatively,  for  small  eccentricities. 


^rp\  Nfr/t,  ^  j  j  b  ~ 
'-p ) 


(c-lU) 


^  1%  iB  used  with  the  eccentric  anwaaly,  E, while  Re  is  used  with  the 
true  anomaly,  ‘U"  ,  After  the  respective  integrations  and  after  small  tenns 
are  dropped,  equal  to  is  found  to  be  consistent  with  the  approxi¬ 
mations. 


1U2 


Figure  C-1. 


Fig.  C-2  shovs  the  results  of  Eq.  (C-lli)  (see  Ref.  l),  as  a  function 
of  the  center  of  mass  position  of  three  elementary  satellite  shapes.  Since 
/  Yp  is  at  most  of  the  order  lO"^,  it.  follows  that  the  parameter 
Arp  rp  Mac  must  be  significantly  less  than  10“1  in  order  for 

the  time  to  dan?)  to  half-amplitude  to  become  a  small  fraction  of  the  decay 
time.  As  seen  from  Fig.  C-2,  the  most  fotrward  locations  of  the  vehicle 
mass  center  produce  damping  times  only  sli^tly  less  than  the  decay  time, 
hence  it  must  be  concluded  that  the  aero<%mamic  danplng  has  little  effect 
on  the  angular  oscillation  over  the  time  period  that  the  orbit  elements 
may  be  regarded  as  constants. 
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